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BEARINGS  ONLY  AIR-TO-AIR  RANGING 


1.  INTRODUCTION 

Airborne  passive  infrared  search  and  track  (IRST)  devices,  once  they  have  detected  and  started 
to  acquire  the  track  of  a  point  target,  do  not  provide  immediate  estimates  of  the  target’s  location, 
range  and  velocity.  It  is  crucial,  of  course,  to  obtain  such  estimates.  As  a  result,  it  is  desirable  to 
have  effective  ranging  algorithms  as  a  component  of  a  passive  IR  detection  system. 

It  is  in  the  nature  of  passive  IR  detection  from  an  ai’bome  platform  that  information  regarding 
die  target  takes  the  form  of  collections  of  bearings  (azimuth  and  altitude  angles)  of  the  target  relative 
to  the  observer,  who  himself  is  in  motion.  Target  parameters,  such  as  location  and  velocity,  are  then 
derived  from  non-linear  Auctions  (typically  as  ratios  of  sums  of  trigonometric  Auctions)  of  the  bear- 


Pasaive  ranging  has  not  yet  received  die  attention  given  to  the  problem  of  target  detection.  As  a 
result,  the  literature  on  the  subject  is  relatively  sparse,  and  the  publications  that  do  exist  [1-3]  deal 
mainly  with  the  simpler  two-dimensional  problem  of  surface-to-surface  passive  ranging.  Furthermore, 
the  focus  of  each  of  the  published  papers  has  been  on  the  presentation  of  a  particular  methodology, 
with  only  a  few  examples  of  its  performance  and  no  real  attempt  to  relate  performance  to  the  parame¬ 
ters  rf  observer  and  target  motion  and  sensor  characteristics. 

The  purpose  of  this  report  is  to  present  the  results  of  a  study  of  the  performance  of  a  variety  of 
air-to-air  passive  ranging  techniques,  particularly  with  regard  to  the  behavior  of  the  mean  and  stan¬ 
dard  deviation  of  range  estimates  as  target  location  and  motion,  observer  motion,  sensor  resolution 
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and  length  of  track  are  varied.  Attention  has  been  concentrated  on  the  cases  of  a  fixed  target  and  of  a 
target  moving  at  constant  velocity,  although  the  computer  programs  developed  to  simulate  and  analyze 
these  cases  can  be  used  with  any  specified  type  of  observer  and  target  motion.  The  methods  selected 
involve  least  squares  fits  using  all  observed  bearings,  least  squares  fits  using  observed  azimuths  only 
(effectively  reducing  the  problem  to  two  dimensions),  and  combining  "minimal  estimates,"  which  are 
closed- form  solutions  for  a  target  that  produce  a  selected  small  number  of  bearings.  Details  of  each 
method  are  presented  in  the  appropriate  sections  of  the  sequel. 

In  this  report,  the  process  of  obtaining  estimates  of  the  range  and,  where  appropriate,  velocity  of 
a  target  by  a  moving  obaerver  it  simulated.  It  is  assumed  that  bearings  are  obtained  at  regular  inter¬ 
vals  (taken  to  be  l  second  in  nearly  all  cases  considered)  and  that  errors  consist  of  independent  Gaus- 
shat  errors  in  and  altitude  measurements  with  zero  mean  and  standard  deviation  of  0.333 

tnrad.  For  fixed  targets,  accurate  estimates  of  target  location  (to  within  a  few  percent  of  the  true 
values)  can  be  obtained  within  about  30  seconds  (assuming  one  set  of  bearings  is  obtained  each 
second)  for  nmrty  any  combination  of  target  location  and  observer  motion.  If  the  target  is  initially 
well  off-axis,  such  accuracy  can  be  obtained  in  considerably  less  tune,  while  if  the  target  is  directly  in 
front  of  the  obaerver  whan  first  detected,  more  time  will  be  needed  for  a  good  estimate.  A  sound 
strategy  then  is  for  die  obaerver,  having  detected  the  target,  to  accelerate  in  order  to  put  the  target  in 
an  off-axis  position  as  quickly  as  possible.  Furthermore,  range  estimates  can  be  obtained  in  a  compu¬ 
tationally  simple  manner  by  using  the  sequence  of  azimuths  to  obtain  successive  estimates  of  the 
target’s  x  and  y  coordinates  and  a  single  altitude  angle  to  provide  an  estimate  of  the  z  component. 
Moving  targets  commonly  require  some  60  seconds  for  good  estimates  of  target  location  and  velocity; 
this  figure  is  quite  sensitive  to  the  specifics  of  the  target’s  initial  position  and  the  relative  motion 
between  observer  and  target.  Unlike  the  fixed  target  case,  where  a  good  strategy  for  the  observer  can 
be  determined  a  priori,  highly  effective  maneuvers  for  the  observer  in  the  case  of  a  moving  target  can 
only  be  determined  after  a  reasonable  estimate  of  the  target’s  velocity  has  been  obtained.  Moreover. 


it  appears  that  least  squares  estimation  of  the  location  of  a  moving  target  will,  for  certain  combina¬ 
tions  of  observer  and  target  motions,  consistently  produce  estimates  based  on  40  to  60  sets  of  bearings 
that  are  for  from  the  true  values.  In  ether  cases,  estimates  cluster,  in  proportions  that  depend  on  tar¬ 
get  and  observer  motion  parameters,  about  two  or  more  widely  separated  values,  one  of  which  may 
be  the  actual  range  of  the  target.  Further  study  is  required  for  a  more  complete  understanding  of  this 


It  should  be  pointed  out  that  the  most  commonly  cited  method  for  performing  passive  ranging 
has  been  the  application  of  an  extended  Kalman  filter  to  the  sequence  of  obtained  bearings.  While 
this  method  has  the  advantage  of  producing  estimates  recursively,  and  therefore  (perhaps)  more 
rapidly,  it  has  the  disadvantages  of  requiring  linearization  of  the  problem,  which  could  entail  a  loss  of 
accuracy  in  die  estimator,  and  of  requiring  initial  estimates  (actually  no  more  than  guesses)  of  target 
parameters  and  estimation  error  statistics,  which  makes  estimates  depend  on  other  quantities  than  the 
observed  bearings.  See  [3]  for  more  on  this  aspect  of  the  Kalman  filter.  Instead,  non-linear  least- 
squares  estimation,  which  in  the  linear  case  is  equivalent  to  Kalman  filtering,  is  used  to  obtain  the 
desired  estimates  in  this  report. 

2.  BACKGROUND 

A  point  target,  located  at  T  -  (x,y,z),  has  bearings  0  and  relative  to  an  observer  at  the  ori¬ 
gin,  where 

tan0-y/x,  sin  $  •  z/r,  r2  -  x2  +  y2  +  z2. 

Bearings-cnly  ranging  methodology  requires  a  moving  observer,  an  appropriate  assumption  for  an 
air-to-air  scenario.  The  coordinate  system  used  will  be  chosen  so  that  at  the  outset  of  the  ranging 
procedure,  when  the  target  is  first  detected  (i.e.,  at  tune  t  »  0),  the  observer  is  considered  to  have  an 
initial  velocity  v0  in  the  +x  direction. 
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In  bearings-only  ranging  the  observer  obtains  a  sequence  of  sets  of  measurements  of  the  target’s 
hearings,  and  from  this  sequence,  together  with  knowledge  of  his  own  trajectory,  attempts  to  deter¬ 
mine  the  relevant  target  motion  parameters,  the  most  important  of  which  are  the  target’s  initial  loca¬ 
tion,  range  and  velocity.  A  single  set  of  bearings  is  not  sufficient  to  estimate  any  of  these  features  of 
the  target’s  motion,  but  it  will  provide  information  about  the  line  of  sight  between  the  observer  and 
the  target  at  the  time  the  bearings  were  obtained.  Additional  sets  of  bearings  then  allow  estimation  of 
target  motion  parameters  through  what  are  effectively  triangulation  procedures;  i.e.,  for  a  fixed  tar¬ 
get,  intersections  of  lines  of  sight  are  determined. 

This  report  will  be  concerned  with  observed  (measured)  bearings,  which  are  estimates  (cor¬ 
rupted  values)  of  the  target's  true  bearings.  Estimates  are,  of  course,  realizations  of  random  variables 
referred  to  as  estimators.  To  be  consistent  with  statistical  usage,  estimators  will  be  denoted  by  upper 
caw  symbols,  whereas  their  estimates  will  be  written  using  the  corresponding  lower  case  characters. 
For  example  estimators  of  8,  the  target’s  azimuth,  and  d,  its  elevation,  may  be  denoted  by  H  and  P 
respectively;  corresponding  estimates  will  be  designated  h  and  p.  In  general,  true  values  of  quantities 
derived  from  9  end  d  will  be  denoted  by  lower  case  Greek  letters,  estimators  by  upper  case  Roman 
letters,  and  measured  values  by  corresponding  lower  case  Roman  symbols.  Frequently,  bearings  will 
be  identified  by  subscripts  or  by  functional  dependence  (0(f),  d(0).  In  such  cases,  H,  P,  and  other 
ttfimtenn  and  all  estimates  will  be  similarly  modified.  If  the  true  bearings  are  6  and  d,  then  the 
measured  bearings,  h  and  p,  aie  assumed  to  be  realizations  of  independent  Gaussian  random  variables 
H  and  P  with  a  common  variance;  i.e.  H  —  N  (0,  <r)  and  P  -  N  (0,  a2).  The  purpose  of  this  report 
is  to  specify  the  statistical  behavior  of  various  bearings-only  ranging  techniques  in  terms  of  observer 
motion,  target  motion  and  measurement  errors. 


4 


3.  LINK  OF  SIGHT  STATISTICS 


Although  a  single  line  off  sight  does  not  specify  anything  about  the  range  of  the  target,  it  does 
provide  an  estimate  of  the  direction  in  which  the  target  lies.  As  a  result,  given  that  the  estimate  of 
the  direction  is  sufficiently  accurate,  all  that  is  needed  to  locate  the  target  is  an  estimate  of  any  of  its 
coordinates.  Prom  the  point  of  view  of  determining  statistical  behavior,  the  process  of  estimating  a 
single  coordinate  of  the  target’s  location  is  much  more  amenable  to  analysis  and  simple  computational 
procedures  than  is  that  of  directly  estimating  its  range.  Present  day  infrared  sensors  provide  estimates 
of  lines  of  sight  that  an  extremely  good,  so  that  we  may  chcoee  to  pursue  bearings-oniy  ranging  as  a 
problem  in  estimation  of  a  single  target  location  coordinate  with  negligible  loss  of  accuracy.  This  is 
illustrated  in  Figure  1,  where  for  a  fixed  "cone”  of  pnobebility,  say  p,  a  single  line  of  sight  is  seen 
to  yield  a  small  spread  compared  to  the  area  common  to  two  such  cones,  which  is  a  region  that  has 
probability  p2  of  containing  the  target  Figure  1  also  illustrates  dun  the  intersection  of  two  lines  of 
sight  produces  a  biased  estimate  of  the  target's  true  position,  and  that  the  estimation  region  becomes 
mote  elongated  aa  the  target's  initial  azimuth  approacha  0  degrees.  Furthermore,  while  the  present 
topic  is  referred  to  as  "passive  ranging,"  it  is  really  "passive  location,"  since  to  know  only  that  a 
potential  target  is,  say,  33  nautical  miles  away  is  to  know  practically  nothing  about  the  position  of  the 
target.  Fortunately,  when  line-of-sigbt  errors  are  small,  the  issues  of  range  and  location  estimation 
are  effectively  equivalent. 

Given  the  true  bearings  9  and  d  of  a  target  relative  to  an  observer  at  the  origin,  the  tme  line  of 
sight  from  the  observer  to  the  target  is  along  the  unit  vector  <  cos 9  cos  d.  sin  9  cos  <£,  sin  <p  > , 
which  we  will  denote  by  <fi,  fj,  $3  > .  On  the  other  hand  observed  bearings  h  and  p  are  realiza¬ 
tions  of  Gatmian  random  variables  H  and  P,  with  H  -  N  (9,  <r)  and  P  ~  N  (9.  <r),  and  they 
determine  an  observed  line  of  sight  with  unit  vector 

<«t,  e3  >  -  <  cos  h  cosp,  sin  h  cos  p,  sinp  >  , 
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which  Sa  in  tun  •  realization  of  the  raodom  vector 


<£|,  £j,  £j>  -  <  cos  H  cos  P ,  sin  H  cos  P ,  sin  P  >  . 

Furthermore,  we  may  write  H  -  9  +  U,  and  P  -  $  +  K,  where  both  U  and  V  are  independent 

Oenaaian  random  variables  with  mean  0  and  variance  a2.  Thus, 

£(  ■  cos  H  coe  P  -  cos  (#  +  £/)  cos  (d  +  10 

*  (cos  9  cos  U  -  sin  •  sin(/)  (cos  d  cos  K  -  sin  <9  sin  K). 

In  practice  we  expect  to  have  a  <  <  1,  so  that  we  may  use  the  approximations  sin  U  =  U  and 

cos  U  *  1  -  U2/ 2,  and  startin':  y  for  V.  Thus 

£  t  *  C(1  -  l/a/2)  cosd-(/sin5][l-  Ka/2)  cos  9  -  Ksin  dl* 

From  the  independence  of  U  end  V  the  expected  value  of  Ex  is  easily  obtained: 

£(£,)  *  cos  I  cos  d  0  -  f/Tp  =*  £,  (1  -  o2). 

After  a  bit  of  algebra  we  also  obtain,  ignoring  terms  higher  than  a2. 

Vat  (£,)  -  £(£?)  -  (£<£,))*  -  (cos2  •  sin2  d  +  sin2  9  cos2  d). 

Sfanflariy, 


£2  ■  sin  cos P  *  1(1  —  (f*/2)  sta  #  +  (/cos  #1  [(l  -  F*/2)cos  d  -  Ksin  dl, 
from  which  it  follows  tbit 

£(J?2>  a  sin# coed  (1  -  a* ft?  *  fed  -  o2),  end 
VariJEj)  »  e2(cos2#cos2d  +  sin2 0 sin3 . 

Although,  £)  can  always  be  obtained  from  the  condition  that  £?  +  £|  +  £|  «  1,  we  still  specify 
the  marginal  distribution  cf  £3  «  sin/*  =  (1  -  Ka/2)sind  +  Kcosd: 

£(£3)  «  siadd  ~  «*/2)  -  fc(l  -  J/T),  and 
K«r(£j)  *  ^co^d- 
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Note  that  tor  certain  values,  e.g.  0  *■  4  «  0,  or  0  «■  •  e/2,  the  shove  formulas  yield  vari¬ 

ances  of  0  in  some  cases.  This  is  of  course  due  to  the  omission  of  high  order  terms.  In  reality,  for 
0  ■  4  m  0  or  w/2,  Kar(£|)  «  a4,  for  example. 

To  obtain  a  more  complete  description  of  the  behavior  of  lines  of  sight,  the  pairwise  behavior  of 
Eif  E2  and  £3  needs  to  be  examined.  Recall  that  the  covariance  of  two  random  variables  X  and  Y, 
denoted  Cov  ( X ,  T),  is  obtained  by  Cov  (X,  Y)  -  E  (XY)  -  E(X)E(Y).  Calculations  yield: 

Cov(Ei£ 2)  *  o*coe0sinFO  —  2cos*d)» 

Cov(Ex,Ey)  =  -  o2 cos 0 cos 4 sin 4,  and 
Cov  (£2,£j)  a  -  a2 sin 0 cos 4 sin 4. 

The  point  of  the  above  analysis  is  that  errors  in  range  estimates  due  to  deviations  of  estimated 
lines  of  sight  from  their  true  values  are  on  die  order  of  a2/?,  where  £  is  the  true  range  of  the  target. 
For  example,  using  a  ■  0.333  mrad,  then  even  for  h  »  0  +  3  <r  and  p  ■»  4  +  3a,  the  range  and 
individual  coordinate  estimates  ate  within  0.1%  of  their  true  values.  Thus  with  negligible  additional 
lots  of  accuracy,  range  estimates  may  often  be  reduced  to  estimates  of  a  single  target  location  coordi¬ 
nate.  Furthermore,  with  knowledge  of  the  individual  and  joint  statistical  behavior  of  the  line  of  sight 
lAMuponents,  one  could  obtain  slight  improvements  on  estimates  based  on  lines  of  sight.  (There 
appears  to  be  little  practical  advantage  in  doing  so,  however.) 

4.  FIXED  TARGET  RANGING 

In  a  number  of  situations  it  is  necessary  to  estimate  the  distance  to  a  target  which  is  known  (or 
assumed)  to  be  stationary.  Examples  are  distant  small  clouds,  nearby  cloud  edges,  or  objects  fixed  on 
the  ground. 


Three  distinct  methods  will  be  presented  for  estimating  the  range  of  a  fixed  target.  The  reason 
for  using  t  multiplicity  of  approaches  is  that  bearings  only  range  estimation  is  a  non-linear  estimation 
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problem*  and  so  the  investigator  is  faced  with  certain  trade-offs  when  selecting  the  animator  to  be 
used  in  a  particular  application.  These  trade-offs  involve  bias  and  efficiency  of  the  estimator,  compu¬ 
tational  requirements  and  ease  of  use  in  recursive  updating  schemes,  and  analytic  tractability.  (In 
linear  estimation  problems,  such  considerations  do  not  arise  since  least  squares  estimators  are  "the" 
choice  according  to  all  of  the  preceding  criteria.) 

All  of  the  methods  operate  on  a  common  set  of  assumptions.  Specifically,  N  sets  of  bearings 
an  assumed,  resulting  in  2N  measurements  (N  tii moths  and  N  altitudes)  with  which  to  estimate  the 
three  target  coordinates.  The  kth  set  of  bearings  in  the  collection  is  denoted  by  h(k)  and  p(k)  and 
the  location  (assumed  known  without  error)  from  which  the  observer  obtained  the  kth  set  of  bearings 
by  (*(*).  y(k),  z(k)).  When  H  is  relevant  to  the  discussion,  we  will  make  use  of  the  fact  that  the 
bearings  are  obtained  in  a  sequential  manner. 

In  assessing  the  performance  of  an  estimator,  many  different  combinations  of  target  location, 
observer  motion  and  other  parameters  of  the  observation  process  must  be  measured.  This  has  been 
done  for  each  of  the  selected  estimators;  results  are  given  in  the  following  sections. 

4.1  Leeet  Sqnares  Range  Estimation 


When  Af  >  1,  the  moat  commonly  used  technique  is  that  of  least  squares  estimation  of  the  unk¬ 
nown  parameters.  For  the  bearings  only  ranging  problem,  two  different  Auctions  suggest  themselves 
as  criterion  Auctions  whose  sum  of  squares  is  to  be  minimized.  The  first  is  simply  the  difference 
between  obeerved  and  fitted  bearings;  i.e.  we  minimize  with  respect  to  x ,  y  and  z: 


S  x 


z 

k 


h(k)  -  Arctan 


[y  -y(*)l 

> 

2  r 

.  T- 

x  -  x (k) 

N.  J 

+  L,  i 

*  L 

p(k)  -  Arcs  in 


r  -\  -\ 

:  -  s(*) 
r(k) 
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where  r2(k)  -  (jc  -  x(k))2  +  (y  —  y(k))2  +  (z  —  z(k))2.  The  second  function  transforms  this 
difference  between  angles  into  a  difference  between  functions  of  the  angles.  Specifically,  we  minim¬ 
ize  with  respect  to  x,  y  and  z : 

*  -  ?  (“*<*>  -  ] : 1  +  5  -  *7^']  ’  ® 

Simulations  have  shown  that  in  many  cases  these  functions  yield  very  similar  estimates,  but  in  those 
cases  where  there  is  a  practical  difference,  minimizing  more  often  produces  estimates  that  are 
closer  to  the  true  target  coordinates.  Consequently,  studies  of  simultaneous  least  squares  estimation 
of  all  three  target  coordinates  have  been  based  upon  minimizing  S\.  A  closed  form  solution  of  the 
normal  equations  for  5t  (i.e.  the  system  of  three  equations  resulting  from  setting  the  partial  deriva¬ 
tives  of  St  with  respect  to  each  of  x,  y  and  z  equal  to  zero)  does  not  appear  to  be  readily  obtainable, 
so  non-linear  least  squares  estimation  software  (the  SNLSE  non-linear  least  squares  estimation  pack¬ 
age  of  the  MATHLJB  subrouting  library)  was  used  to  produce  estimates  from  simulated  sample  data. 

In  order  to  assess  the  performance  of  this  least  squares  estimator,  a  variety  of  initial  target  loca¬ 
tions  and  observer  motions  of  either  constant  velocity  or  constant  acceleration  were  selected.  For 
each  such  target/observer  combination  the  number  of  bearings  taken  and  the  time  between  successive 
bearings  were  varied  methodically.  With  all  relevant  parameters  chosen,  simulations  were  performed 
by  determining  the  correct  bearings,  corrupting  tnem  by  addition  of  independent  samples  from  a 
Gaussian  N  (0,  a2)  distribution,  and  then  using  the  non-linear  least  squares  software  to  obtain  an  esti¬ 
mate  of  the  target  location.  Typically  this  was  repeated  100  times,  so  that  for  a  selected  scenario  the 
statistics  (mean,  median,  standard  deviation,  third  moment,  and  selected  percentiles  for  each  target 
coordinate  as  well  as  range)  for  a  sample  of  size  100  were  determined.  Tables  A1  and  A2  (see 
Appendix)  show  the  results  of  two  such  analyses,  while  Fig.  2  is  a  plot  of  the  sample  standard  devia¬ 
tions  as  a  function  of  the  number  of  points  on  the  track,  for  a  selected  set  of  values  of  time  between 
bearings.  Examination  of  the  tables  and  the  plot  supports  the  intuitively  appealing  notion  that  for  a 
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Nusbsr  of  Points  on  Track 


given  amount  of  time,  the  estimator  tha»  involves  the  largest  number  of  bearings  provides  the  best 
estimate  of  target  location.  Consequently,  all  further  analyses  were  carried  out  with  bearings  being 
obtained  at  the  rate  of  one  set  per  time  unit  (assumed  to  be  1  second).  Further  examination  of  Tables 
A1  and  A2  also  supports  the  argument  based  on  the  geometry  of  Figure  1  that  range  estimates  are 
biased  towards  overestimating  the  true  range. 

Even  after  the  interval  between  successive  bearings  has  been  fixed,  a  large  number  of  combina¬ 
tions  of  sensor,  observer  and  target  parameters  still  remain.  In  order  to  reduce  the  number  of  cases 
to  a  manageable  one,  while  preserving  the  geometric  flavor  of  the  relation  between  observer  motion 
and  target  position,  two  particular  quantities  were  kept  fixed  for  most  of  the  analyses.  These  quanti¬ 
ties  are  the  standard  deviation,  a,  of  the  measurement  errors  in  the  target’s  azimuth  and  altitude, 
which  was  held  at  0.333  mrad,  and  the  initial  observer  velocity,  which  was  fixed  at  C.16  nm/s  in  the 
x  direction.  These  values  were  chosen  because  0.333  mrad  is  a  realistic  value  for  errors  in  modem 
passive  infrared  sensors,  while  0.16  nm/s  represents  a  typical  speed  for  a  patrol  aircraft.  Two  exam- 
plea  of  the  effects  of  varying  the  sensor's  resolution  are  given  in  Figures  3  and  4.  These  plots  and 
their  accompanying  tables  (Tables  A3  and  A4  in  the  Appendix)  imply  that,  with  all  other  parameters 
held  fixed,  the  standard  deviation  of  the  least  squares  range  estimator  varies  approximately  linearly 
with  the  standard  deviation  of  sensor  error. 

The  plots  that  follow  will  be  of  standard  deviations  of  range  estimates  (or  their  logarithms) 
versus  number  of  points  on  the  track,  with  the  latter  being  equivalent  to  the  length  of  time  the  target 
is  being  tracked.  Standard  deviations  have  been  chosen  because  they  are  the  most  common  measure 
by  which  to  express  expected  variability.  Recall  that  for  a  Gaussian  random  variable  approximately 
68%,  93%  and  99.3%  of  sample  values  will  fall  within  1,  2,  and  3  standard  deviations,  respectively, 
of  the  mean  of  the  random  variable.  Although  least  squares  range  estimators  do  not  appear  to  be 
Gaussian,  particularly  for  short  tracks,  for  sufficiently  long  tracks  (of  length  at  least  20)  the  above 
percentages  appear  to  be  reasonable  estimates.  For  example,  based  on  Table  A2  the  statement  can  be 
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made  wish  a  high  degree  of  assurance  that,  at  least  90%  of  the  range  estimates  of  a  target  initially 
located  at  the  point  (50,  40,  2)  run  would  be  within  0.5  tun  of  the  true  range  of  64.06  run  for  an 
observer  traveling  with  initial  velocity  <.16,  0,  0>  nm/s  and  constant  acceleration  <0,  -.006,  0> 
nm/s2,  obtaining  bearings  every  second  for  25  seconds.  Also,  although  least  squares  estimators  are 
biased,  die  bias  has  been  found  to  be  small  compared  to  expected  sampling  error,  so  that  nearly  all  of 
the  information  about  die  reliability  of  least  squares  estimators  lies  in  the  sample  standard  deviation. 

With  the  resolution  of  the  sensor  and  the  initial  velocity  of  the  observer  fixed,  the  remaining 
relevant  parameters  are  the  target's  range,  its  initial  direction  relative  to  the  observer,  and  the  motion 
of  the  observer.  As  for  as  initial  target  direction  and  observer  motion  are  concerned,  two  characteris¬ 
tic  cases  will  be  considered  for  each:  (a)  targets  that  are  initially  “in  front  of*  the  observer;  i.e.  tar¬ 
gets  which  have  initial  azimuth  of  approximately  10  degrees  or  less,  and  (2)  targets  that  are  “well 
off-axis'’  —  targets  with  initial  azimuth  approximately  40  degrees  or  more.  It  is  tacitly  assumed  that 
the  sensor  is  limited  to  azimuths  -45  and  45  degrees,  although  on  occasion  this  limitation  is  ignored. 
The  observer  was  taken  to  be  moving  with  either  constant  velocity  or  constant  acceleration.  In  the 
latter  case  accelerations  of  approximately  lg  (.006  nm/sec2)  were  used  in  the  simulations.  For  targets 
initially  in  front  of  the  observer,  two  different  constant  accelerations  were  considered,  while  for  off- 
axis  targets  only  a  constant  velocity  observer  and  one  with  constant  acceleration  in  the  positive  x 
direction  are  treated.  The  reason  for  the  latter  restriction  is  that  once  the  target  is  well  off  axis,  any 
appreciable  acceleration  in  any  but  the  x  direction  will  either  remove  the  target  from  the  sensor’s  field 
of  view  or  bring  the  target  into  a  less  favorable  head-on  position. 

The  effect  of  range  alone  is  seen  in  Figures  5  through  9,  where  for  fixed  observer  motion  and 
target  direction  we  see  the  effects  of  ever  increasing  range.  Figures  5  and  6  deal  with  a  target  ini¬ 
tially  well  off-axis  (initial  azimuth  ■  38.7  degrees)  for  a  constant  velocity  and  accelerating  observer 
respectively,  while  Figs.  7,  8  and  9  repeat  the  situation  for  a  head-on  target  (initial  azimuth  *  1 1.3 
degrees).  Note  that  these  are  plots  of  the  logarithms  of  sample  standard  deviations  against  track 
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Imgth.  Logarithm  plots  won  chosen  in  order  to  (It  •  wide  range  of  values  on  a  single  plot,  but 
incidentally  the  plots  show  that  from  about  20  seconds  on  there  is  a  roughly  log-linear  relation 
between  track  length  and  sample  standard  deviation.  This  suggests  that  after  approximately  20 
secoode,  estimation  error  falls  off  approximately  exponentially  with  length  of  track.  Examination  of 
tha  tabiaa  tor  those  graphs  (Tables  A7  through  All)  indicates  that,  for  any  of  the  plots,  doubiini  the 
target  range  roughly  multiplies  the  sample  standard  deviation  by  four.  This  suggests  that  the  standard 
deviation  of  tha  learn  squares  range  eerimmor  varies  as  tha  square  of  the  range,  independent  of  track 

The  preceding  flee  graphs  were  founded  to  show  the  effects  of  range  on  accuracy  of  the  estima¬ 
tor.  The  next  two  graphs  fodfoais  tha  effect  of  imkial  target  direction  on  the  estimate.  Figures  10  and 
11  (else  Tabiaa  A 12,  A13)  presses  the  errors  in  animating  the  range  of  a  target  initially  40  nm  away 
Anna  tha  oberrer  for  bridal  asknutha  of  43*.  40*.  30*.  20*,  13\  10*  and  3*  for  two  different 
observer  modona.  Norn  that  far  an  aooaforedag  observer,  tha  effect  of  initial  target  direction  essen¬ 
tially  dhappama  after  23  aecxmdi,  Figure  12  (afoo  Tibiae  A14,  A15)  shows  the  results  of  varying  the 
dhecdoe  of  a  target  whose  initial  x  oowdlimhi  is  fixed,  in  this  instance  at  30  nm.  Note  that  for  short 
trader,  target  direcdoe  affects  dnminaW;  i.e.  for  up  to  13  seconds  earimetaa  of  the  target  at  (30.  10, 
2)  (Rage  •  31  am)  are  lire  relfobia  than  those  of  the  target  at  (30,  40,  2)  (Range  -  64  nm). 
dmpilB  the  feet  font  hi  the  tetter  case  tha  target  is  13  am  Anther  away.  For  sufficiently  long  tracks, 
hammer,  range  pfays  the  doednmtt  rale.  Figures  13  and  14  are  dose-up  looks  at  the  curves  in  Figure 
12. 


The  fled  p— mm  is  observer  motion.  Figures  13  through  17  (where  Fig.  16  is  t  doseup  of 
tha  right  end  of  Pig.  15)  show  dm  effects  of  varying  observer  acceleration  for  on-  and  off-axis  tar¬ 
gets.  Norn  that  in  tha  short  ran  very  unrelfobic  eadmams  can  arise  from  an  observer  accelerating  in 
tha  "wrong”  direction,  while  over  a  longer  period  of  time  acceleration  in  any  direction  produces 
more  aocurase  range  animates  than  does  remaining  at  the  same  velocity.  This  of  course  is  the  result 
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of  the  observer  moving  so  that  the  target  passes  directly  in  front  of  him  somewhere  along  his  trajec¬ 
tory.  Not  included  here  are  plots  of  other  possible  observer  motions,  such  as  acceleration  out  of  the 
x  -  y  plana.  Such  accelerations  were  in  fact  considered,  but  all  cases  examined  strongly  suggest  that 
the  only  important  features  of  observer  acceleration  are  the  magnitude  of  the  acceleration  in  the  x 
direction  and  the  magnitude  orthogonal  to  the  x  direction. 

In  summary,  the  salient  features  of  least  squares  estimation  of  position  and  range  of  a  stationary 
target  are:  (1)  the  most  reliable  mtimams  come  from  targets  that  are  well  off  axis,  (2)  for  short 
tracks  reliability  is  itasamihmri  more  by  the  direction  of  the  target  than  by  its  range,  while  the  situa¬ 
tion  reverses  itself  for  longer  tracking  times,  (3)  regardless  of  target  direction  and  range,  after  some 
30  seconds  the  standard  deviation  of  the  range  estimate  is  no  mote  than  a  percent  or  two  of  the  true 
range  of  the  target. 

While  the  least  squares  estimator  presented  here  offers  reliable  estimates  within  a  reasonably 
short  period  of  dan,  the  method  has  some  distinct  drawbacks.  The  first  disadvantage  is  that  of 
repairing  computationally  complex  non-linear  least  squares  software.  While  such  software  is  readily 
available,  its  use  consumes  considerable  amounts  of  compmer  time.  For  example,  it  required  an 
average  of  one  and  a  half  hours  of  CPU  time  on  a  VAX  11/780  to  yield  the  168  sample  standard 
deviations  that  were  plotted  in  each  of  Figures  3  through  9.  Another  disadvantage  of  this  approach  is 
that  it  does  not  take  advantage  of  the  sequential  manner  in  which  bearings  are  obtained  by  the 
observer.  An  ideal  scheme  would  use  each  new  pair  of  bearings  to  update  the  current  estimate  of  tar¬ 
get  location  in  a  computationally  simple  manner.  (A  truly  ideal  method  would  be  simple  enough  to 
allow  updating  to  be  done  by  hand,  particularly  during  test  and  debugging  phases.)  For  linear  estima¬ 
tion  problems  the  desired  adaptive  estimation  procedure  is  available  in  the  form  of  the  Kalman  filter, 
which  is  a  recursive  method  for  determining  the  standard  linear  least  squares  fit  of  the  model  to  the 
data.  Kalman  filters  have  been  applied  to  ship- to- ship  passive  ranging  problems  [2,3],  where  the 
problem  is  viewed  as  two-dimensional,  bearings  consist  only  of  azimuths,  and  the  model  has  been 
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linearized  in  order  to  apply  the  filter.  One  of  tfaeae  references  [3]  dwells  on  the  pitfalls  of  this 
method,  particularly  retarding  the  sensitivity  of  short  track  estimates  to  the  values  of  required  initial 
estimates  that  must  be  supplied  before  any  measurements  have  actually  been  taken.  Fundamental  to 
the  use  of  the  Kalman  filter  is  an  initial  —«*"*■»*  of  the  error  covariance  matrix  for  the  quantities 
being  estimated,  in  addition  to  an  initial  of  their  values.  Unfortunately,  the  covariance 

depends  on  the  target  s  range,  so  initial  estimates  are  apt  to  be  very  unreliable,  and  it  may  take  quite 
a  long  time  for  the  data  to  correct  for  a  poor  guess.  On  the  other  hand,  non- recursive  least  squares 
fitting  requires  only  an  initial  position  estimate,  to  which  the  procedure  is  not  particularly  sensitive. 
A  final  drawback  of  both  is  that  they  are  purely  numerical  schemes,  providing  no  analytic  insight  into 
the  roiss  played  by  target  range  and  direction  or  by  observer  motion  in  the  statistical  behavior  of  the 


Since  it  has  already  bean  eantbHahad  that,  given  small  errors  in  bearings  measurements,  accurate 
osrimathm  of  a  single  target  cnorriinata  is  wfflrtent  to  yield  good  estimates  of  the  other  coordinates, 
and  since  much  of  the  difficulty  in  dealing  analytically  with  the  and  S2  criterion  functions  arises 
from  terms  involving  altitude  angles,  azimuth-only  estimates  of  x  and  y  were  considered.  These 
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However,  Sj  does  not  lend  itself  to  a  simple  closed  form  solution  of  its 
does  S4,  but  a  modified  form  of  it  does  provide  the  desired  simplicity. 


(3) 

(4) 

normal  equations.  Neither 
Setting  g(k)  -  tan(A(*)), 


define 


Ss  -  El(y  -  y(*))  -  g(k)  (x  -  x(k » i2. 

k 


(5) 
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««ff— *«t  to  zero  the  partial  derivatives  of  Ss  with  respect  10  x  and  y ,  we  obtain  the  normal  equations 
for  Sy 

y  £f  <*>  “  £* (*)y<*)  -  x  £**<*)  +  £**(*)»(*)  -  o 


Ny  -  £y(*)  -  X Et(k)  +  £/(*)»(*)  -  0. 

Solving  these  equations,  we  obtain aotatione Jt  and?  that  satisfy 

yZt(k)  -  f  £#(*)»(*)  -  Zt <*fr(*)  +  tiHk)x(k) 
Lt\k)  -  i£g(k) 


(6) 


$  -y  +  4t  - (7) 

where  jr  -  ^£y(Jt) and#  -  ^ £#(*). 

These  estimate,  which  oaa  be  eery  easily  updated  as  new  aahanths  see  obtained,  appear  deceptively 
stasflar  an  eetimases  of  slope  and  intercept  la  Uaser  tegreaaioa  problems.  Here  however,  unlike  the 
Unoar  regrseeioa  case,  i  and  f  are  non  Hauer  fonctiona  of  the  observations,  which  we  may  take  to  be 
the  g  (k)’s,  end  so  their  statistical  behavior  cannot  be  readily  dstenninsd  in  aa  analytic  maimer. 
Instead,  as  was  done  In  the  previous  section,  shnulations  are  performed.  Once  again  sensor  resolution 
is  set  at  0.333  mrad,  initial  observer  velocity  at  0.16  nm/s  in  the  x  direction,  and  the  same  combina¬ 
tions  of  target  range  and  direction  and  observer  motion  are  uaad.  In  addition,  the  same  sequence  of 
rfWit"—  "mt*—  (using  IMSL’s  GGNML  generator  of  values  from  the  standard  normal  distribution) 
is  employed  in  performing  both  sets  of  simulations.  This  azimuth-only  procedure  yields  estimates  of 
jt  end  y  only.  Fnimarton  of  target  range  requires  an  estimate  of  the  target’s  z  -coordinate  as  well. 
Bated  on  the  argument  drat  a  single  line  of  sighs  is  e  very  good  estimator  of  the  true  target  direction, 
aatimtmr  of  z  were  obtained  via  the  formula 

i  -  (tanp(O))^ (i2  +  y2),  (8) 
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where  />(0)  is  tbs  dmtoi  observed  at  time  0.  The  results  of  these  simuletioQS  are  presented  in 
Table  A16.  For  comparison,  the  corresponding  results  using  the  standard  least  squares  approach  of 
the  previous  section  can  be  found  in  Table  A7.  From  these  tables  it  is  clear  that  «zimutb-only  least 
squares  estimation  tranches  the  accuracy  of  estimation  using  all  obtained  bearings.  Also  of  interest  is 
the  feet  that  for  short  tracks  azimuth-only  estimates  tend  to  underestimate  the  target's  true  range.  An 
***** fanrina  of  the  procedure  with  o  set  to  aero  shows  that,  for  a  typical  ranging  scenario,  the  equa¬ 
tions  for  i  sad  S  do  indeed  yield  x  sad  y  ea  solutions,  but  they  are  not  especially  well  conditioned. 
Uhl  suggests  that  the  “»*■— -ti— in»t  of  range  might  be  a  computational  artifact,  so  the  procedure 
was  npsaesd  using  double  precision  arithmetic.  This  produced  mean  ranges  that,  surprisingly,  were 
even  forther  from  the  true  range,  aftek  by  rathur  small  amounts  and  with  the  anticipated  smaller  stan¬ 
dard  tr*— Tbs  insrrnhfehla  Implication,  la  that  the  oberved  underestimation  is  the  result  of  bias 
in  tha  tathuatnr.  Further  study  of  the  hies  in  azimuth-only  estimators  is  proceeding.  Azimuth-only 
eedamthm  has  bean  briefly  mmdoned  hi  the  Uaerature  [5],  where  k  is  compared  with  the  Kalman 
fllMr  hi  fen  feipto-eMp  ranging  pmhlmn  (whim  feet*  an  only  azimuths)  and  where  k  appears  in  e 
rather  different  fliinnilwina  It  is  of  significant  inaareet  that  the  avenge  emonnt  of  CPU  time  required 
to  gaosnua  fee  dm  for  ae  azimuth-oaty  eathnrinn  equivalent  to  that  in  Figures  5  through  9  was 
•haul  seven  nriuens,  compered  wife  fee  hour  end  e  half  needed  for  the  least  squares  method  using  ell 


In  summary,  azimuth-only  rang*  estimation  is  compondonally  very  efficient  and  practically 
nuerhns  the  accuracy  of  the  toon  complicated  standard  knit  squares  method.  Its  main  drawbacks  are 
due  k  appears  to  email  a  bin  of  pntwfly  unknown  nature  and  once  again  is  not  sufficiently  amenable 
to  analysis  to  make  k  ezpUck  how  the  mean  cud  variance  of  the  range  estimator  depend  on  target 
location,  obeerver  motion  and  sensor  resolution,  hi  addition,  azimuth-based  methods  do  not  apply  in 
the  general  moving  target  cam. 
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Band  oa  MUmI  Estimators 


The  third  method  for  estimating  range  presented  in  this  report  aims  at  uncovering  the  manner  in 
which  target  and  observer  parameters  affect  the  statistics  of  the  estimator,  at  the  expense  of  increased 
uncertainty  in  the  estimator.  It  is  intended  not  so  much  as  a  practical  estimation  scheme  (at  least  in 
the  fixed  target  case,  where  the  azimuth-only  least  squaies  approach  is  preferable),  but  as  a  means  of 
throwing  light  or  he  nature  of  the  biases  and  variability  inherent  in  bearings  only  estimation.  The 
basic  idee  is  to  use  only  as  many  of  the  bearings  as  are  required  in  order  to  obtain  a  dosed  form  solu¬ 
tion  for  the  target  location  that  yielded  those  bearings.  Because  the  resulting  range  estimator  uses  the 
possible  number  of  bearings,  it  is  referred  to  ns  »  “minimal**  estimator.  Since  minimal 
range  estimators  turn  out  to  have  simple  formulations  in  terms  of  the  bearings  that  define  them,  it  is 
possible  to  obtain  formulas  for  their  means  and  variances  that  hold  to  very  good  approximations. 
Being  baaed  on  er,  few  bearings,  minimal  estimators  are  of  course  not  very  accurate.  Accuracy  can 
be  increased,  and  use  made  of  all  obtained  sets  of  bearings,  by  combining  collections  of  minimal  esti¬ 
mates  into  a  single  better  ettimaSe. 

In  the  statin— ty  target  case  the  three  target  coordinates  are  to  be  estimated,  so  that,  at  the 
minimum,  two  sets  of  bearings  must  be  obtained.  Two  seta  of  bearings  yield  four  measurements,  and 
since  they  are  corrupted  values  of  the  true  bearings,  they  produce  an  ovetdetermined  system.  There 
is  a  choice  of  methods  for  dealing  with  this  situation.  Having  four  measurements  available  to  esti¬ 
mate  three  location  coordinates  allows,  by  selecting  all  possible  sets  of  three  measurements,  four  dis¬ 
tinct  minimal  estimates  of  the  target  location.  These  estimates  may  be  considered  singly,  or  as  a 
composite  «■**»«—»  such  eg  a  or  median. 

While  the  two  sett  of  bearings  to  be  taken  could  be  obtained  from  any  two  points  in  space. 
intended  applications  of  this  study  as  well  as  avoidance  of  computational  complexity  suggest  that  the 
fint  bearings,  (Ao,Pq),  be  taken  from  the  origin,  arid  the  second  set,  (hy,px),  from  the  point  (d,0,0). 
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This  choice  of  Mcood  location  is  also  justified  by  the  feet  that  minimal  estimates  are  often  quick  esti¬ 
mates,  and  u  such  are  besed  on  bearings  taken  in  rapid  succession.  As  a  result,  since  we  are  assum¬ 
ing  the  x-axis  of  our  coordinate  system  to  be  oriented  with  the  observer's  initial  velocity,  we  expect 
that  in  a  short  time  interval  only  a  net  change  in  the  observer's  x  coordinate  can  realistically  be 
achieved. 


With  the  observer  motion  specified  in  the  previous  paragraph,  and  with  realistic  target  locations 
relative  to  an  awumed  airborne  observer,  various  estimation  procedures  were  examined  for  the  pur¬ 
pose  of  selecting  one  for  use  in  the  ensuing  analysis.  While  the  performances  of  the  various  methods 
were  expected  to  be  comparable  for  a  target  well  off-axis,  it  is  desired  that  the  chosen  estimator  per¬ 
form  notably  better,  in  terms  of  both  bias  and  variance,  when  the  off-axis  azimuth  6®  is  near  zero. 
Conceptually,  the  chosen  estimation  procedure  uses  A®  and  A|  to  obtain  an  azimuth-only  estimate  of 
the  target’s  x  and  y  morrfinntes,  and  than  uaet  p®  to  obtain  a  farther  estimate  of  the  z  coordinate, 
whaaoe  the  range  be  estimated. 

As  meatkxted  at  the  beginning  of  this  report,  the  target  is  initially  located  at  (x,y,z)  relative*  to 
the  observer.  Hearing!  A®,  p®,  and  At,  with  A|  being  the  target’s  azimuth  relative  to  (<f  ,0,0),  yield  a 
unique  nartmnte  of  the  target  location,  t  ■  (f ,  obtained  as  follows.  Letting  o®  »  tan  A®  and 

«i  •  tan  A(  we  may  writs 


and 

Solving  fori  and?  yields 


and 


a®  ■  tan  A®  -  j>/i 
o,  -  tan  A!  -  ^/(i  -  d), 

i  -  da{/(ax  -  <*o) 


9  ■  a®?  »  dflofl,/(a ,  -  ao). 


(9) 

(10) 

(U) 

(12) 
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Using  these  satlmarea  and  the  line  of  sight  determined  by  Hq  and  pg  we  can  determine  !.  Letting 

h0  -  sec  A0tanp0  (13) 

it  follows  from  projecting  (£,£,£)  into  the  x-y  plane  that 

I  “  bof  m  dh0fl|/(al  —  flo).  (14) 

Thus  the  estimated  target  position  is  7  -  (i,?,£)  -  i<  l,a0,h0>.  Now,  as  discussed  in  the  sec- 

don  on  the  statistics  of  Unas  of  sight,  for  the  small  values  of  <r  under  consideration,  <  l,a0,b0>  is  a 

very  good  tsthnsis  of  <  1,  taaflo,  sec  9otan4o> .  the  true  vector  relating  the  coordinates  of  the  target 

relative  to  the  origin.  Errors  in  estimating  7,  then,  arise  mainly  from  errors  in  £.  Consequently, 

attention  will  be  concentrated  on  the  analysis  of  the  random  variable  £,  particularly  with  regard  to  its 

depandance  on  both  T  and  d. 

Note:  The  preceding  discussion  does  not  apply  directly  when  x  *  0,  for  then  both  a0  and  b0 
might  be  arbitrarily  large.  In  this  case,  we  would  simply  solve  the  same  equations  explicitly  for  9, 
and  obtain  7  -  $  <a0\  1,  h0'>,  for  appropriate  a0*  and  b0'»  “d  perform  a  similar  analysis  for 
This,  however,  will  not  be  e  problem  here,  since  we  are  assuming  that  |  0q|  does  not  exceed  45 
degrees. 

&  is  given  by  X  -<£4,/Mi  -  Atf,  where  A9  -  tanH0>  «nd  Ax  -  tanH,.  Now 
Ho  -  +  V.  end  H,  -  d|  +  K,  where  is  the  true  azimuth  of  the  target  from  the  origin,  0]  is 
the  true  azimuth  from  (4,0,0),  and  U  and  V  are  independent  Gaussian  variables  with  mean  0  and  a 
common  variance  o2(o2<  <  1).  Now  tantf0  -  tan(0o  +  U)  (tan0o  +  tan£/)/(l  -  tmtitfanU).  Let¬ 
ting 


ocq  -  tan0o  ■  y/x  .  (15) 

and  recalling  that  for  small  tan  ^  we  have 


<40  ~  (ao  +  £/)/(!  -  o of/). 


(16) 
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Maflariy,  defining 


a,  ■  tan  0,  ■  y/{x  -  d). 


Ax  *  (o,  +  V)/(\  -  a, 10. 


(17) 


(18) 


Then 


Thb 


i  - 


rf(al  +  IQ 
1  -  g|K 


(4,  -  a,  +  V  a*  +  V 

1  -  a,K  “  1  -  acotf 


(19) 


d{*x  +  K)  (1  -  d0U) 


(at  -  a*)  +  (1  +  citfax)  (V  -  U)  +  (a,  -  oto)  UV 
Factoring  (a,  -  <%j)  .«t  of  tbe  denominator,  and  setting 


(20) 


K  - 


«i  -  ca# 


C  - 


1  +  «o«i 
<*\  -  <*0 


we  obtain 


(21) 


JT(a,  +  V)  (1  -  «9lO 


(22) 


Tba  following 
of  t,  aod  are  immediate  from  the  deftaitioaa: 


1  +  C(K  -  £0  +  UV ' 
among  C,  K,  oq,  alf  jc,  y  and  d  will  be  uaeftil  in  determining  the  behavior 


k  .  m.-i) , 

£«i  “  x,  /Too  •  x  -  d  , 

c  *(x  -  <0  +  K+JL 

yd  d  ■ 
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To  obtain  estimates  of  E(X)  and  Var(X)  that  are  accurate  up  to  the  <?  term,  note  that  since  U 
and  V  are  independent  and  each  has  mean  zero,  the  lowest  order  term  to  which  UV  can  make  a  con¬ 
tribution  is  the  o4  term.  Thus  we  may  ignore  UV  terms,  and  write 


a  _  £(«i  ~  otpai  U  +  V ) 

X  *  1  +  C(V  -  U) 

Under  the  approximation  1/(1  +  w)  a  1  -  w  +  w2,  this  becomes 


(24) 


X  a  K(ax  -  <*<pt\U  +  U)(  1  -  C(V  -  U)  +  C\V  -  P)2). 

Performing  the  specified  multiplication,  and  ignoring  terms  whose  expected  values  are  of  degree  three 

and  higher  in  <r,  we  obtain 

£  a  K(ctx  -  otxC(V  -  U)  +  axC\V  -  U)2  -  cl&lxU  +  a<fitxCU(V  -  U)  +  V  -  CV(V  -  U )) 
Taking  expected  values  and  recalling  that  E((V  —  U)2)  »  2a2, 


EQt)  «  K(ax  +  -  «o a,  -  1)) 

*  Kax  +  a2C(2CKax  —  Kct(fitx  —  AT).  (25) 

Since  Kax  »  x  and  y  »  oto*,  this  becomes 


E(X  ~  x  +  azC(2xC  -  y  -  K). 

Using  the  relation  K  +  y  *  Cd,  we  have 

£(X)  a  x  +  t^C2^  -  d).  (26) 

Squaring  £,  dropping  higher  order  and  mixed  terms,  taking  expected  values  and  using  the  relations 
(23)  once  again  we  obtain 

E(X2)  a  x2  +  +  y1  +  K2  -  4 a0Cx2  -  4 KCx). 

Then 


Var{X)  -  E(X2)  -  E(X)2  ~ 

fQCh2  +  lC2xd  +  y2  +  K2  -  Ac^Cx2  -  4 KCx).  (27) 
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Using  the  relations  K  »  Cd  -  y  and  y2  *  Cdy  -  x(x  -  d),  this  simplifies  to 

ttrrtf)  a  ^[C^lc2  -  2xd  +  d2)  -  2x{x  -  d)].  (28) 

Remark:  In  most  applications  x>>d,  so  with  very  little  additional  loss  in  accuracy  the 
moments  of  £  take  on  much  simpler  forms;  viz. 

E(X)  *  x(l  +  2<j2C2)  (29) 

Vartf )  *  or2[2r2(C2  -  1)].  (30) 

Furthermore,  in  this  case,  C2>  >  1,  so  that  we  may  use 

Kirtf)  *  2x2CV.  (31) 

Note  that  for  a  fixed  range  and  displacement  d,  C  increases  as  y  decreases,  i.e.  as  the  target 

becomes  more  “head-on.”  As  a  result  both  the  bias  and  the  variance  of  X  increase.  Not  only  that, 
but  the  second  order  approximation  used  for  1/(1  4-  C(V  -  (/))  becomes  less  accurate  as  C 
increases.  Thus  even  the  large  values  of  EQt)  and  V ar(jt)  obtained  above  are  actually  underesti¬ 
mates  of  the  true  values  when  C  is  large.  That  C  becomes  infinite  when  the  target  is  directly  in  front 
of  the  observer  is,  of  course,  a  consequence  of  the  assumption  that  the  second  set  of  bearings  is  also 
taken  from  a  point  on  the  x-axis.  In  practice,  one  need  not  take  the  second  bearings  on  the  x-axis, 
but  in  realistic  situations  the  point  from  which  they  are  taken  will  be  quite  near  the  x-axis.  As  a 
result,  there  will  be  increases  in  both  the  bias  and  variance  of  AT  as  an  estimator  of  x  when  range 
increases,  when  d  decreases,  and  for  fixed  range  and  d,  when  y  decreases. 

To  improve  estimates  of  E(jt)  and  Var(X)  when  C  is  large,  fourth  order  approximations  of  X 
were  obtained.  Using  1/(1  +  w)  »  1  —  w  +  w2  -  w3  +  w4,  including  UV  terms,  performing  the 
required  algebra,  noting  that  those  terms  in  which  C  has  the  same  exponent  as  a  dominate  the  others, 
taking  expected  values  and  again  assuming  that  x  >  >d  and  C2  >  >  1 ,  the  results  are 


E(X)  *  x(l  +  2o2C2  +  \2a*C*) 

(32) 

Var(X)  ~  2x2CV  +  32*2CV. 

(33) 
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The  preceding  estimates  are  the  results  of  a  great  many  truncations,  and  it  is  reasonable  at  this 
point  to  inquire  as  to  how  good  they  are  in  practice.  Table  1  shows  the  estimated  means  and  standard 
deviations  using  both  the  second-order  and  fourth-order  approximations  to  X  for  a  few  cases.  Also 
included  in  this  table  are  the  results  of  simulations  in  which  the  means  and  standard  deviations  of 
samples  of  size  1000  were  obtained.  Note  that  for  large  values  of  C  even  the  fourth  order  approxi¬ 
mation  may  grossly  underestimate  the  sample  standard  deviation.  This  arises  from  the  fact  that  in 
reality  X  has  infinite  mean,  and  that  in  any  finite  sample  soi  e  extremely  large  or  even  negative 
values  are  likely  to  occur.  Fortunately  these  extreme  values  do  not  correspond  to  realistic  situations, 
so  in  practice  such  an  estimate  can  be  ignored  if  it  arises.  With  this  in  mind,  extreme  scores  were 
effectively  eliminated  by  excluding  the  upper  and  lower  1%  of  the  scores  in  the  samples  with  large 
C.  The  statistics  of  these  middle  98%  of  the  samples  are  shown  in  parentheses  in  the  table,  where 
appropriate.  With  this  proviso  we  see  that  for  targets  well  off  axis  or  with  sufficient  separation 
between  bearings,  the  second  order  approximations  are  good  ones,  while  in  other  cases  including  the 
fourth  order  term  produces  better  estimates. 


Table  1  —  Minimal  (Two-point)  Estimation  of  Fixed  Target  Position 
Bearings  taken  from  (0,  0,  0)  and  (d,  0,  0) 

Sample  values  are  from  samples  with  N  -  1000  and  a  »  0.333  mrad 
Sample  values  in  parentheses  refer  to  middle  98%  of  scores 
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In  practice,  of  course,  one  does  not  expect  to  be  restricted  to  two  sets  of  bearings  in  estimating 
the  location  of  the  target.  Instead,  a  collection  of  sets  of  bearings  would  be  available,  the  collection 
being  obtained  in  a  sequential  manner,  with  new  pairs  of  bearings  being  added  at  regular  intervals.  A 
method  was  sought  to  select  pairs  of  sets  of  bearings  and  combine  the  resulting  minimal  estimates  into 
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a  single  effective  range  estimate,  preferably  a  scheme  that  can  be  easily  updated  as  new  sets  of  bear¬ 
ings  are  obtained.  Examination  of  sampling  distributions  of  minimal  estimators  indicate  that  the  sam¬ 
ple  median  is  a  much  better  estimator  of  the  true  target  location  parameters  than  is  the  sample  mean, 
especially  for  small  separations  between  the  points  from  which  the  pair  of  bearings  is  taken.  (That 
the  median  is  a  better  estimator  than  the  mean  also  follows  from  the  fact  that  the  form  of  the  estima¬ 
tor  of  x  is,  effectively,  x/(l  +  W)  where  W  is  normally  distributed  with  mean  0.  Since,  on  the 
average,  half  of  the  W’s  in  a  sample  are  positive,  and  the  other  half  negative,  half  of  the  estimates 
will  be  less  than  x ,  and  half  will  be  larger.  Thus  the  sample  median  is  an  unbiased  estimator  of  x , 
while  the  mean  of  x/(l  +  W)  is  approximately  x(l  +  Var{W)).)  The  procedure  chosen  takes  a 
sequence  of  sets  of  bearings,  forms  minimal  estimates  of  the  target’s  x  coordinate  from  each  pair  of 
azimuths  using  formula  (6),  and  then  takes  the  median  of  this  collection  of  minimal  estimates  to  be 
the  estimate  of  x .  Estimates  of  y,  z,  and  the  range  are  then  obtained  using  this  x  estimate  and  the 
initial  line  of  sight  This  procedure  is  still  under  investigation,  particularly  regarding  the  best  choice 
of  pairs  of  points  to  be  used,  so  results  are  still  preliminary.  Table  A17  gives  the  results  of  using  this 
method  in  two  familiar  cases. 

5.  CONSTANT  VELOCITY  TARGETS 

Conceptually,  the  methods  for  estimating  die  target  motion  parameters  of  a  moving  target  are 
the  same  as  those  for  a  fixed  target.  The  practical  difference  is  that  there  are  more  parameters  to  be 
estimated,  just  how  many  mote  depending  on  the  target  motion  model  selected.  For  instance,  if  a 
constant  velocity  target  is  assumed,  there  are  three  velocity  components  to  be  estimated  in  addition  to 
three  location  coordinates,  while  if  a  target  with  constant  acceleration  is  assumed,  there  are  yet 
another  three  acceleration  components  to  be  estimated.  Specifying  that  the  target  move  along  a  partic¬ 
ular  type  of  parametric  path,  e.g.  along  the  arc  of  a  circle,  requires  estimating  as  many  quantities  as 
are  needed  to  completely  specify  the  motion.  The  effects  of  having  more  parameters  to  estimate  are 
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that  more  complicated  relations  exist  between  observer  and  target  motions  and  the  reliability  of  esti- 
mators,  that  longer  tracks  must  be  obtained  in  order  to  estimate  range  and  other  quantities  to  specified 
accuracies,  and  that  the  vagaries  of  non-linear  estimation  procedures  with  their  reliance  on  numerical 
methods  must  be  accepted. 

In  this  section  attention  will  be  restricted  to  constant  velocity  targets.  Only  one  estimation 
method  will  be  examined  in  detail— least  squares  fitting  of  the  constant  velocity  model  to  the  observed 
bearings.  In  this  case,  there  are  six  quantities  to  estimate,  the  three  target  coordinates,  x ,  y  and  z ,  as 
well  as  its  velocity  components  denoted  by  Vx ,  Vy  and  Vz . 

Studies  of  alternative  ranging  methods  are  ongoing.  Unfortunately,  moving  targets  do  not 
necessarily  stay  at  the  same  altitude,  so  that  the  azimuth-only  least  squares  method  that  works  so  well 
for  fixed  targets  is  not  applicable  to  moving  targets  in  general.  However,  in  many  cases  a  distant  tar¬ 
get  will  not  undergo  significant  changes  in  its  altitude  during  the  period  of  its  observation,  so  that 
azimuth-based  estimation  procedures  are  currently  being  considered.  Also,  as  in  the  fixed  target  case, 
minimal  closed-form  estimates  of  target  motion  parameters  are  available,  in  this  case  requiring  three 
sets  of  bearings  to  provide  sufficient  data  to  estimate  the  six  target  values.  Collections  of  such  esti¬ 
mates  may  be  combined  to  form  a  single  improved  value.  Concurrent  with  present  efforts  in  the 
fixsd  target  case,  the  emphasis  of  ongoing  work  in  minimal  estimation  of  moving  target  parameters  is 
on  the  optimal  choice  of  sets  of  bearings  to  be  combined  for  a  single  minimal  estimate. 


5.1  Least  Squares  Range  Estimation 

A  similar  procedure  is  used  as  in  the  fixed  target  case,  based  on  the  differences  between 
observed  and  fitted  bearings;  the  expression 
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it  minimized  with  respect  to  x,  y,z.  Vx,  Vy,  Vz, 


where  r*(k)  -  (x  +  (k  -  l)Kc  -  x(jfc))2  +  (y+  (k  -  1  )Vy  -  y(k))2  +  (z  +  (k  -  1)K? 

-  zto)2. 

Once  again  the  properties  of  the  estimator  were  studied  by  examining  its  performance  under  a 
wide  variety  of  simulated  scenarios,  and  again  the  standard  deviation  of  the  100  range  estimates 
obtained  for  each  chosen  combination  of  observer  and  target  parameters  was  used  as  the  primary 
measure  of  reliability.  In  ranging  a  constant  velocity  target,  the  observer  cannot  himself  travel  at 
constant  velocity,  for  then,  even  in  the  absence  of  any  measurement  errors,  the  obtained  sequence  of 
bearings  does  not  yield  a  unique  solution  for  the  target  perimeters.  This  can  be  very  easily  seen  by 
"Mgwwf  two  targets,  each  flying  parallel  to  the  observer,  the  first  target  at  half  the  range  and  speed 
of  the  second.  As  long  as  the  obaerver  maintains  the  same  velocity,  these  targets  will  yield  the  same 
sets  of  Consequently,  in  all  cases  ***«■■"*»  we  have  assumed  an  accelerating  obaerver, 

■orf,  farther,  that  the  acceleration  is  W  **—*•"» 

Figures  IS  through  21  show  the  standard  deviation  of  range  estimates  as  a  function  of  track 
length  for  a  variety  of  combinations  of  observer  and  target  motions  (also  see  Tables  A18  and  A19). 
In  these  figures  trade  lengths  begin  at  15  (still  assuming  one  set  of  bearings  per  second)  since  for 
shorisr  tracks  the  eethneses  are  so  poor  that  they  provide  no  usefal  information.  Figure  18  deals  with 
a  target  that  is  initially  in  front  of  the  obaerver,  while  Figs.  19,  20  and  21  examine  an  off-axis  target. 
Even  a  cursory  examination  of  these  graphs  reveals  a  great  deal  more  complexity  when  the  target's 
velocity  aunt  be  nitimaswd  in  addition  to  ha  range.  The  most  apparent  difference  is  that  the  standard 
deviation  of  the  range  aatimawu  docs  not  monotonically  decrease  with  the  length  of  the  track.  Indeed 
there  are  often  dramatic  rises  in  the  sample  standard  deviation.  In  an  effort  to  understand  this 
phenomenon,  percentiles  were  obtained  for  range  estimates  from  samples  of  size  1000  simulating  a 
target  originally  at  (50,  40,  1)  nm  with  velocity  <>.20,  -  20,  0)  run/s  being  tracked  by  an  observer 
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Nuabsr  of  Points  on  Track 


Least  Squares  fit  to  a  Constant  Velocity  Target 
St.  Dev.  of  Range  vs  *pts  on  track 
0  samples/point,  Initial  Observer  Vo  -  <.16,0,0> 
Target  at  (50,40,1)  —  R  -  64.04 


r* 
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Number  of  Points  on  Track 
Figure  19 
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Munbcr  of  Point*  on  Track 

Figure  21 


with  acceleration  <0,  *.006,  .006>  nm/s3.  Ten  such  percentile  distributions  were  constructed,  and 
in  each  it  was  found  that  there  appear  to  be  two  distinct  clusters  of  estimates,  with  approximately 
15%  of  the  range  estimates  being  near  49  nm  and  the  remainder  scattered  within  a  few  nautical  miles 
of  65  nm.  The  sample  data  are  listed  in  Table  A20. 

Figures  18  through  21  are  the  usual  graphs  of  sample  standard  deviations  versus  track  lengths, 
and  while  they  inti  irate  irregularity  in  the  estimation  process,  they  do  not  explain  its  source.  Figure 
22,  which  is  typical  of  what  can  arise  in  non-linear  least  squares  estimation,  provides  some  insight 
into  the  nature  of  the  irregularities  in  figs.  18  thorough  21.  In  Figure  22  the  mean  of  the  100  esti- 
maass  is  plotted  as  a  function  of  track  length,  together  with  error  bars  corresponding  to  one  sample 
tfaadud  deviation.  Here  the  mean  range  drifts  away  from  the  true  range  for  a  rather  long  interval  of 
track  lengths.  Thus  a  bias,  which  can  be  much  larger  than  the  sample  standard  deviation  and  which 
dtp—**  on  track  length  m  well  as  target  and  observer  motion  parameters,  is  introduced  when  target 
velocity  moat  alao  be  mtimaaad.  As  in  Table  A20,  there  may  be  two  (or  more)  well  separated  values 
around  which  efimmi  cluster,  this  accounts  for  the  unusually  large  sample  standard  deviations  found 
for  tang  tracks.  This  phenomenon  is  most  likely  due  to  the  feet  that  in  a  non-linear  estimation  prob¬ 
lem  m  six  dimensions,  rather  different  combinations  of  target  location  and  velocity  may  yield  very 
similar  sequences  of  bearings  relative  to  the  specified  observer  trajectory.  Measurement  errors  then 
remit  in  random  selection,  according  to  some  as  yet  undetermined  rules,  of  one  of  the  possible 


Further  analysis  is  being  performed  to  gain  a  more  thorough  understanding  of  the  combined 
effects  of  target  location  and  velocity  and  observer  motion  on  the  performance  of  least  squares  estima¬ 
tors,  both  in  the  case  of  a  constant  velocity  target  and,  in  what  promises  to  be  much  more  compli¬ 
cated,  the  situation  of  an  accelerating  target.  A  final  graph,  Figure  23,  compares  the  results  when  a 
fixed  target  is  known  to  be  fixed,  and  so  can  be  ranged  using  the  methods  of  section  4,  versus  when 
it  is  not  assumed  to  be  stationary,  and  so  its  lack  of  motion  must  be  inferred  from  its  bearings.  This 
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Nunber  of  Points  on  Track 

Figure  22 
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of  Point*  on  Track 


wggerts  that  a  comparative  study  of  the  ranging  of  slowly  moving  targets  by  both  methods  might  be 
ftuMhl. 


6.  SUMMARY  AND  CONCLUSIONS 

The  petfbnnanoe  of  various  air-to-air  passive  ranging  techniques  has  been  examined  as  a  func¬ 
tion  of  target  location  and  motion,  observer  motion,  and  length  of  track.  Bearings  are  obtained  at 
ona  second  intervals  and  there  are  independent  Gaussian  errors  in  azimuth  and  altitude. 

For  fixed  targets,  three  range  estimation  methods  have  been  considered,  each  generally  giving 
accurate  estimates  of  target  location  for  track  lengths  of  30  seconds  or  less: 

1)  Least  Squares  Range  Estimation— This  method  provides  accurate  range  estimates  within  a 
short  tracking  time,  but  requires  computationally  complex  software  and  extensive  computer  processing 
time;  it  doee  not  update  target  location  estimates  using  sequential  bearings. 

2)  Azimuth  Only  Leatt  Squares  Range  Estimation— This  method  is  nearly  as  accurate  as  the  pre¬ 
vious  method  bet  is  computationally  much  mote  efficient.  Target  location  estimates  are  readily 
updated  as  new  azimuths  are  obtained.  This  method  is  the  best  of  the  three  methods  for  the  fixed  tar¬ 
get  case  but  is  not  directly  adaptable  to  the  moving  target  case. 

3)  Range  Estimation  Baaed  on  Minimal  Estimators— This  method  uses  only  a  small  number  of 
bearings  so  tint  a  dosed  form  solution  for  the  target  location  can  be  found.  Consequently,  range  esti¬ 
mates  are  not  as  accurate  as  in  die  previous  methods,  but  explicit  information  can  be  obtained  on  the 
effect  of  target  and  observer  parameters  on  the  estimator  statistics. 

For  moving  targets,  a  track  length  of  60  seconds  is  typically  required  for  good  estimates  of  tar¬ 
get  location  and  velocity,  and  this  figure  is  highly  dependent  on  the  observer/target  geometry.  Least 
Squares  Range  Estimation  has  been  investigated  for  the  constant  velocity  target  case.  However,  it 
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was  found  that  in  some  cases  range  estimates  have  a  cluster  point  separate  from  that  at  the  true  range. 
Further  study  will  be  required  to  develop  a  technique  to  resolve  this  ambiguity. 
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APPENDIX 


Table  A1 


Least  Squares  Estimates  -of  Fixed  Target  Position 
100  Samples  per  Cell  a  -  0.333  mrad 

Target  Location:  (50,50,2)  Range  -  70.7*1  no 


Constant 

Observer  Velocity:  <.16,0,0> 

nm/s 

Time  Between 

Bearings  (,*«*> 

1 

2 

4 

6 

8 

#  of  Bearings 

2 

83.63 

73.09 

70.94 

71  .84 

71 .21 

61  .80 

12.65 

5.23 

3.62 

2.71 

3 

73.26 

70.68 

70.66 

70.85 

70.74 

12.31 

5.41 

2.45 

1  .73 

1  .23 

4 

71  .67 

70.86 

70.98 

70.69 

70.62 

6.76 

3.23 

1  .57 

1.19 

0.83 

5 

70.81 

70.76 

70.93 

70.76 

70.69 

4.75 

2.23 

1  .23 

0.88 

0.54 

6 

70.52 

70.70 

70.70 

70.74 

70.71 

3.45 

1  .65 

0.81 

0.57 

0.41 

8 

70.87 

70.74 

70.66 

70.70 

70.73 

1  .96 

0.97 

0.55 

0.32 

0.27 

9 

70.71 

70.77 

70.71 

70.77 

1 .85 

0.90 

0.42 

0.32 

10 

70.99 

70.79 

70.75 

70.77 

1  .59 

0.71 

0.39 

0.24 

12 

70.73 

70.79 

70.74 

1.17 

0.59 

0.26 

15 

70.67 

70.73 

70.71 

0,96 

0.34 

0.20 

16 

70.88 

70.75 

0.73 

0.44 

18 

70.72 

70.76 

0.60 

0.31 

20 

70.69 

70.73 

0.52 

0.26 

24 

70.73 

0.37 

25 

70.71 

-  0.40 

30 

70.74 

0.31 

The  upper 

entry  In 

each  cell,  is  the  sample  mean, 

while  the 

lower 

entry  is 

the  sample 

standard  deviation. 
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Table  A2 

Least  Squares  Estimates  of  Fixed  Target  Position 

|  100  Samples 

per  Cell 

o  -  0.333 

or  ad 

Target  Location:  (50,40,2) 

Range  -  64. 

06  nm 

Initial 

Observer  Velocity: 

<  .  1 6 , 0 , 0>  nm/s 

Observer  Acceleration:  <0,- 

.006, 0>  nm/s 

2 

Time  Between  Bear  ings  (&£) 

1 

2 

4 

6 

8 

#  of  Bearings 

2 

75.85 

66.11 

64.53 

64.98 

64.42 

62.07 

9.71 

4.91 

2.93 

2.11 

3 

66.32 

64.00 

63.98 

64.14 

64.06 

10.88 

4.58 

1  .88 

1  .24 

0.83 

4 

64.87 

64.14 

64.23 

64.02 

63.99 

5.89 

2.63 

1.13 

O'.  77 

0.49  i 

5 

64.08 

64.07 

64.19 

64.08 

64.04  j 

• 

4.02 

1  .73 

0.83 

0.53 

0.30 

6 

63.90 

64.03 

64.05 

64.06 

64.04 

2.86 

1  .23 

0.51 

0.31 

0.12 

8 

64.15 

64.06 

64.01 

64.04 

64.06 

1  .55 

0.68 

0.31 

0.16 

0.12 

9 

64.04 

64.07 

64.05 

64.08 

1.44 

0.60 

0.23 

0.15 

10 

64.25 

64.09 

64.06 

64.08 

1 .21 

0.46 

0.20 

0.11 

12 

64.06 

64.09 

64.07 

0.86 

0.36 

0.13 

15 

64.01 

64.06 

64.05 

0.67 

0.19 

0.09 

16 

64.16 

64.07 

0.50 

0.24 

18 

64.05 

64.08 

0.39 

0.16 

20 

64.03 

64.06 

0.34 

0.13 

24 

64.06 

0.23 

25 

64.04 

0.24 

30 

64.06 

0.17 

The  upper  entry  In 

each  oell  is  the  sample  mean,  while  the 

lower 

•ntry  la 

! 

the  sample 

standard  deviation 

e 
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TABLE  A3 


Least  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Cell  Various  o 
Target  Location:  (50,40,2)  Range  -  64.06  nm 

Initial  Observer  Velocity:  <.16,0,0>  nm/s 
Observer  Acceleration:  <Q,-.006,0>  nm/s* 


#  Points 

o  -  .  1 

mr  ad 

o  -  . ; 

333  -irad 

o  ■  .  I 

566  mrad 

on  Traok 

X 

s 

X 

3 

X 

s 

4 

64.27 

1  .77 

64.18 

5.58 

67.97 

16.78 

5 

64.23 

1  .22 

63.74 

4.03 

66.19 

9.18 

6 

64.00 

0.95 

65.28 

2.93 

64.19 

6.54 

8 

64.05 

0.54 

63.81 

1  . 80 

64.14 

3.73 

9 

64.10 

0.43 

64.12 

1.39 

64.43 

2.92 

10 

64.13 

0.37 

64.10 

1  .20 

64.55 

2.53 

12 

64.03 

0.27 

64.04 

0.86 

63.87 

1  .79 

15 

64.06 

0.19 

63.95 

0.60 

64.07 

1  .24 

16 

64.03 

0.16 

63.95 

0.54 

63.87 

1  .04 

18 

64.07 

0.13 

64.11 

0.49 

64.10 

0.86 

20 

64.03 

0.11 

64.05 

0.33 

63.84 

0.74 

24 

64.05 

0.07 

64.04 

0.27 

64.02 

0.50 

25 

64.06 

0.07 

64.07 

0.23 

64.02 

0.47 

30 

64.05 

0.05 

64.05 

0.15 

63.99 

0.35 

32 

64.06 

0.04 

64.06 

0.15 

64.08 

0.25 

35 

64.06 

0.03 

64.07 

0.12 

64.05 

0.23 

36 

64.06 

0.03 

64.06 

0.1  1 

64.07 

0.23 

40 

64.06 

0.03 

64.08 

0.08 

64.07 

0.20 

45 

64.06 

a.  02 

64.06 

0.08 

64.07 

0.13 

48 

64.06 

0.02 

64.06 

0.07. 

64.07 

0.12 

50 

64.06 

0.02 

64.06 

0.06 

64.08 

0.12 

55 


Table  A4 

Least  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Cell  Various  o 
Target  Location:  (50,10,2)  Range  -  51.03  nm 

Constant  Observer  Velocity:  <.16,0,0>  nm/s 


•'  • 

o  -  .1 

mrad 

o  ■  . 

333  mrad 

o  ■  . 

666  mrad 

#  Points 
on  Track 

X 

s 

X 

s 

X 

s 

4 

51  .65 

3.50 

55.57 

16.03 

186.7 

799.7 

5 

51  .18 

2.56 

52.53 

8.86 

57.88 

23.03 

6 

51.12 

1  .85 

51  .90 

6.99 

56.42 

29.67 

8 

50.99 

1  .28 

51  .1  1 

4.27 

51  .97 

8.97 

9 

51  .00 

1  .08 

51.14 

3.69 

51  .82 

7.93 

10 

51  .07 

0.70 

51  .14 

2.40 

51  .68 

5.07 

12 

51.10 

0.57 

51.30 

1  .93 

51  .74 

4.00 

15 

50.97 

0.46 

50.87 

1  .52 

50.82 

3.08 

16 

51 .1  2 

0.41 

51  .35 

1  .39 

51  .74 

2.83 

18 

51  .00 

0.35 

50.96 

1.16 

50.94 

2.33 

20 

51  .03 

0.28 

51  .05 

0.95 

51.11 

1  .90 

24 

51  .07 

0.22 

51.17 

0.72 

51  .33 

1  .44 

25 

51  .03 

0.19 

51  .05 

0.65 

51  .08 

1 .30 

30 

51  .04 

0.16 

51  .08 

0.52 

51.15 

1  .04 

32 

51  .01 

0.13 

50.98 

0.44 

50.93 

0.89 

35 

51  .03 

0.12 

51  .04 

0.41 

51  .05 

0.81 

36 

51  .02 

0.13 

51  .00 

0.42 

50.98 

0.85 

40 

51  .03 

0.09 

51  .04 

0.29 

51  .05 

0.58 

45 

51  .04 

0.08 

51  .07 

0.25 

51.12 

0.50 

48 

51  .03 

0.06 

51  .04 

0.20 

51  .05 

0.39 

50 

51  .03 

0.07 

51  .02 

« 

0.22 

51  .01 

0.44 
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Table  A5 


Least  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Cell  Various  o 
Target  Location:  (50,40,2)  Range  -  64.06  nm 

Constant  Observer  Velocity:  <.16,0,0>  nm/s 


Points 

c  *  .1 

arad 

o  ■  .  ] 

333  mrad 

o  -  .  i 

666  mrad 

Track 

X 

s 

X 

3 

X 

s 

4 

63.84 

1  .70 

64.74 

6.31 

64.48 

12.48 

5 

64.20 

1  .40 

63.11 

4.18 

66.17 

9.94 

6 

64.18 

0.92 

63.95 

2.83 

65.42 

6.62 

8 

64.08 

0.61 

64.09 

2.28 

64.38 

4.06 

9 

64.08 

0.49 

63.91 

1  .65 

64.30 

3.28 

10 

64.00 

0.40 

64.30 

1  .50 

63.70 

2.65 

1  2 

64.05 

0.30 

63.86 

1  .01 

64.04 

2.00 

15 

64.07 

0.24 

64.16 

0.79 

64.17 

1  .60 

16 

64.1  1 

0.22 

64.06 

0.66 

64.42 

1  .46 

18 

64.05 

0.16 

64.02 

0.59 

64.01 

1  .09 

20 

64.05 

0.16 

64.12 

0.54 

64.02 

1  .03 

24 

64.06 

0.13 

64.03 

0.35 

64.08 

0.86 

25 

64.08 

0.10 

64.02 

0.39 

64.17 

0.69 

30 

64.06 

0.08 

64.08 

6.29 

64.05 

0.52 

32 

64.06 

0.06 

64.07 

0.22 

64.05 

0.43 

35 

64.06 

0.06 

64.04 

0.22 

64.08 

0.41 

36 

64.05 

0.06 

64.05 

0.21 

64.00 

0.38 

40 

64.06 

0.05 

64.08 

0.16 

64.04 

0.33 

45 

64.06 

0.04 

64.06 

0.14 

64.03 

0.29 

48 

64.06 

0.04 

64.06 

0.13 

64.05 

0.23 

50 

64.07 

0.03 

64.06 

0.11 

64.09 

0.22 

57 


Table  A6 


Laaat  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Cell  Various  o 
Target  Location*  (50,10,2)  Range  -  51.03  nm 

Initial  Observer  Velocity:  <.16,0,0>  nm/s 
Observer  Acceleration:  <0,-.006,0>  nm/s* 


#  Points 

o  ■  .  1 

nr  ad 

o  ■  . 

333  mrad 

o  -  .  i 

666  mrad 

on  Track 

X 

s 

X 

'  3 

X 

s 

4 

51  .47 

3.02 

51  .27 

10.87 

106.0 

379.1 

5 

50.86 

1  .75 

52.46 

6.11 

52.58 

15.96 

6 

51  .06 

1.34 

51.30 

4.89 

52.77 

10.67 

8 

50.95 

0.71 

51.10 

2.46 

50.91 

4.86 

9 

51  .12 

0.57 

51.13 

2.15 

51  .91 

3.90 

10 

51  .07 

0.48 

50.99 

1  .37 

51  .44 

3.19 

12 

51  .00 

0.33 

51  .15 

1  .08 

50.92 

2.23 

15 

51  .00 

0.21 

51  .08 

0.61 

50.84 

1 .40 

16 

51  .05 

0.18 

51  .06 

0.62 

51  .20 

1.18 

13 

51  .03 

0.13 

51  .04 

0.46 

51  .04 

0.86 

20 

51  .02 

0.10 

51  .02 

0.32 

50.97 

0.68 

24 

51  .03 

0.06 

51  .06 

0.23 

51  .02 

0.43 

25 

51  .02 

0.06 

51  .04 

0.19 

51  .00 

0.42 

30 

51  .03 

0.04 

51  .05 

0.14 

51  .04 

0.29 

32 

51  .03 

0.03 

51  .05 

0.12 

51  .06 

0.22 

35 

51  .03 

0.03 

51  .03 

0.09 

51  .03 

0.20 

36 

51  .03 

0.03 

51  .03 

0.09 

51  .03 

0.17 

40 

51  .03 

0.02 

51  .03 

0.06 

51  .02 

0.14 

45 

51  .03 

0.01 

51  .04 

0.05 

51  .02 

0.09 

48 

51  .03 

0.01 

-51  .03 

0.04 

51  .02 

0.10 

50 

51  .03 

0.01 

51  .03 

0.04 

51  .03 

0.07 

38 


Table  A7 


Least  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Cell  a  .  0.333  mrad 

Constant  Observer  Veloolty*  <.i6,0,0>  nm/s 


T» ( 20 , 

16.2) 

T-(30, 

24.2) 

T  » (  40  , 

32.2) 

T  « ( 50 . 

.40.21 

Points 

R_-  25 

.69 

R_-  38 

1.47 

R_-  51 

.26 

R  -  64.06 

n  Track 

X 

s 

X 

s 

X 

s 

x 

3 

4 

25.83 

0.99 

38.77 

2.10 

51  .40 

3.75 

63.72 

5.76 

5 

25.63 

0.63 

38.55 

1.54 

51  .75 

2.89 

64.76 

4 . 72 

6 

25. TO 

0.52 

38.53 

1.10 

51.18 

2.26 

64.58 

3.15 

8 

25.65 

0.31 

38.44 

0.77 

51  .26 

1.34 

64.15 

2.02 

9 

25.74 

0.26 

38.45 

0.66 

51  .39 

1  .08 

64.14 

1  . 62 

10 

25.71 

0.23 

38.49 

0.43 

51  .43 

0.96 

63.86 

1  .33 

12 

25.68 

0.17 

38.52 

0.34 

51.18 

0.72 

64.04 

1  .00 

15 

25.67 

0.12 

38.44 

0.27 

51  .26 

0.52 

64.11 

0.80 

16 

25.71 

0.11 

38.53 

0.25 

51  .18 

0.44 

64.23 

0 . 73 

18 

25.69 

0.08 

38.45 

0.21 

51  .28 

0.38 

64.03 

0.55 

20 

25.68 

0.07 

38.47 

0.17 

51.16 

0.33 

64.04 

0.52 

24 

25.69 

0.05 

38.50 

0.13 

51  .24 

0.24 

64.07 

0.43 

25 

25.69 

0.05 

38.47 

0.12 

51  .24 

0.23 

64.12 

0 . 35 

30 

25.69 

0.04 

38.48 

0.10 

51.23 

0.18 

64.06 

0.26 

32 

25.70 

0.03 

38.46 

0.08 

51  .28 

0.13 

64.06 

0.21 

35 

25.69 

0.03 

38.47 

0.07 

51  .26 

0.12 

64.07 

0.21 

36 

25.69 

0.03 

38.47 

0.08 

51  .26 

0.12 

64.03 

0.19 

40 

25.69 

0.02 

38.47 

0.05 

51  .27 

0.11 

64.05 

0.17 

45 

25.69 

0.02 

38.48 

0.05 

51  .27 

0.08 

64.05 

0.15 

48 

25.69 

0.02 

38.47 

0.04 

51  .27 

0.07 

64.06 

0.12 

50 

25.69 

0.01 

38.47 

0.04 

51.27 

0.07 

64.08 

0.11 

T"(60, 

48.2) 

T- (80 , 

64.2) 

T-C100 

.80.2) 

T-(120 

.96  2  1 

Points 
n  Track 

R_-  76 
X 

.86 

s 

R  -  102 

X 

.47 

s 

R_-  128.08 
x  s 

153.69 

X  A 

4 

78.00 

9.21 

100.96 

17.23 

128.57 

24.88 

159.52 

41.31 

5 

76.97 

6.04 

103.70 

1 1  .24 

132.02 

17.11 

154.62 

29.96 

6 

76.73 

4.08 

104.03 

8.42 

129.21 

14.71 

163.26 

20.83 

8 

76.93 

3.30 

102.10 

5.26 

128.45 

8.43 

152.62 

12.15 

9 

76.66 

2.38 

103.06 

4.28 

128.60 

7.95 

1 54 . 37 

9 . 72 

1  0 

77.22 

2.18 

102.68 

3.83 

128.04 

*.91 

154.26 

8 . 62 

1’2 

76.58 

1 .46 

102.67 

2.93 

128.70 

4.59 

1 53.63 

6.22 

15 

77.00 

1.15 

102.58 

2.17 

128.30 

2.99 

152.90 

4.62 

1  6 

76.85 

0.95 

102.58 

1.93 

128.22 

3.20 

152.85 

4 .29 

18 

76.81 

0.86 

102.42 

1.43 

128.13 

2.47 

154.14 

4.05 

20 

76.95 

0.78 

102.43 

1.34 

128.06 

1  .96 

153.62 

2.82 

24 

76.81 

0.52 

102.63 

1  .03 

128.31 

1.57 

153.56 

2.42 

25 

76.81 

0.57 

102.45 

0.97 

128.16 

1 . 30 

153-75 

2.18 

30 

76.89 

0.42 

102.41 

0.72 

128.29 

1.12 

153.57 

1  .57 

32 

76.87 

0.32 

102.47 

0.64 

128.23 

0.98 

153.68 

1.48 

35 

76.83 

0.32 

102.57 

0.52 

128.11 

0.91 

153.79 

1  . 27 

36 

76.84 

0.31 

102.48 

0.62 

128.13 

0.85 

153.67 

1.19 

40 

76.88 

0.24 

102.50 

0.42 

1  28.05 

0.61 

153.86 

0.92 

45 

76.87 

0.21 

102.48 

0.39 

128.21 

0,63 

153.71 

0.92 

48 

76.86 

0,19 

102.47 

0.34 

128.14 

0.53 

153.66 

0.81 

50 

76.86 

0.17 

102.49 

0.29 

128.12 

0.52 

153.68 

0.68 

59 


Table  A 8 


I 

i 

i 


Least  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Call  j  -  0.333  arad 

Initial  Observer  Velocity:  C.16,0,0>  no/s 
Observer  Acceleration:  <.006,0,0>  nm/s* 


T-(20, 

16,2) 

T-(30. 

24,2) 

#  Points 

R  -  25 

.69 

R  -  38 

.47 

on  Traok 

X 

s 

X 

3 

25.82 

0.99 

38.75 

1 .98 

5 

25.69 

0.59 

38.59 

1  .44 

6 

25.69 

0.98 

38.52 

1  .00 

8 

25.66 

0.27 

38.99 

0.68 

9 

25.73 

0.23 

38.95 

0.57 

10 

25.71 

0.20 

38.99 

0.36 

12 

25.68 

0.19 

38.51 

0.28 

15 

25.68 

0.10 

38.95 

0.21 

16 

25.70 

0.08 

38.96 

0.19 

18 

25.69 

0.06 

38.97 

0.16 

20 

25.69 

0.05 

38.99 

0.12 

21 

25.69 

0.03 

38.97 

0.09 

25 

25.69 

0.03 

38.48 

0.08 

30 

25.69 

0.02 

38.47 

0.06 

32 

25.69 

0.02 

38.47 

0.05 

35 

25.69 

0.02 

38.47 

0.04 

36 

25.69 

0.01 

38.47 

0.04 

*0 

25.69 

0.01 

38.47 

0.03 

<5 

25.69 

0.01 

38.47 

0.02 

48 

25.69 

0.01 

38.47 

0.02 

50 

25.69 

0.00 

38.47 

0.02 

T-(60, 

98,2) 

T« ( 80 , 

64,2) 

#  Points 

R_-  76 

.86 

R  -  102 

.47 

on  Traak 

X 

a 

X 

4 

4 

77.9a 

8.69 

100.94 

16.41 

5 

76.92 

5.61 

103.51 

10.38 

6 

77.79 

3.73 

103.76 

7.86 

8 

76.91 

2.90 

102.09 

4.62 

9 

76.67 

2.06 

102.97 

3.71 

10 

77.17 

1  .85 

102.62 

3.23 

12 

76.63 

1  .20 

102.63 

2.41 

15 

76.97 

0.90 

102.55 

1.72 

16 

76.85 

0.79 

102.54 

1 .50 

18 

76.82 

0.65 

102.44 

1 .08 

20 

76.93 

0.57 

102.44 

0.99 

29 

76.82 

0.35 

102.57 

0.71 

25 

76.83 

0.38 

102.45 

0.66 

30 

76.88 

0.26 

102.43 

0.45 

32 

76.86 

0.19 

102.46 

0.39 

35 

76.89 

0.18 

102.53 

0.30 

36 

76.85 

0.18 

102.48 

0.36 

90 

76.88 

0.13 

102.49 

0.24 

95 

76.86 

0.10 

102.48 

0.20 

98 

76.86 

0.09 

102.48 

0.17 

50 

76.86 

0.08 

102.47 

0.14 

T-(40, 

32,2) 

T- ( 50 , 

40,2) 

R  -  51 

.26 

R  -  64 

.06 

X 

3 

X 

s 

51  .37 

3.52 

63.72 

5.42 

51  .72 

2.67 

64.66 

4.39 

51.17 

2.06 

64.52 

2.86 

51  .25 

1.18 

64.14 

1  .78 

51  .36 

0.94 

64.13 

1  .42 

51  .40 

0.81 

63-87 

1.14 

51  .19 

0.59 

64.04 

0.82 

51  .26 

0.40 

64.10 

0.63 

51 .20 

0.34 

64.19 

0.56 

51 .27 

0.28 

64.04 

0.41 

51  .19 

0.24 

64.04 

0.37 

51  .25 

0.16 

64.07 

0.29 

51  .25 

0.15 

64.10 

0.23 

51  .24 

0.11 

64.06 

0.16 

51.27 

0.08 

64.06 

0.13 

51  .26 

0.07 

64.07 

0.12 

51  .26 

0.07 

64.04 

0.11 

51 .27 

0.06 

64.06 

0.09 

51 .26 

0.04 

64.06 

0.07 

51  .27 

0.03 

64.06 

0.06 

51.27 

0.03 

64.07 

0.05 

T-(100 

,80,2) 

T-(120 

.96.2) 

R_-  128.08 

R  -  153.69 

X 

s 

X 

s 

128.62 

25.21 

159.01 

38.09 

131 .57 

15.69 

154.26 

26.74 

128.89 

13.20 

162.40 

18.83 

128.35 

7.46 

152.59 

10.66 

128.48 

6.84 

154.24 

8.41 

128.02 

4.26 

154.10 

7.31 

128.55 

3.76 

153.64 

5.18 

128.26 

2.35 

153.04 

3.64 

128.19 

2.48 

153.03 

3.32 

128.12 

1  .87 

154.02 

3.01 

128.05 

1.42 

153.62 

2.06 

128.23 

1  .08 

153.58 

1  .67 

128.13 

0.89 

153.71 

1  .45 

128.21 

0.72 

153.60 

0.97 

128.18 

0.61 

153.68 

0.92 

128.09 

0.53 

153.74 

0.76 

128.11 

0.49 

153.67 

0.70 

128.07 

0.33 

153.78 

0.50 

128.15 

0.32 

153-69 

0.47 

128 .1 1 

0.26 

153.67 

0.41 

128.09 

0.25 

153.68 

0.34 

60 


Tabic  A 9 


Least  Squares  Estimates  of  fixed  Target  Position 
100  Saaples  per  Cell  o  *  0.333  orad 

Constant  Observer  Velocity:  <.16.0,0>  nm/a 


T-(20, 

4.2) 

T-(30, 

6.2) 

T- ( 40 , 

8.2) 

T-( 50 , 

10,2) 

Points 

8-20 

.49 

R  -  30 

.66 

R  -  40 

.84 

R  -  40 

.03 

i  Track 

X 

s 

X 

s 

X 

s 

X 

3 

1 

20.85 

1  .84 

31  .73 

4.34 

41.72 

7.79 

52.43 

14.69 

5 

20.40 

1.14 

30.95 

2.98 

42.29 

5.42 

53.64 

10.34 

6 

20.57 

0.93 

30.82 

2.21 

40.88 

4.60 

52.39 

6.59 

6 

20.45 

0.52 

30.64 

1.46 

40.95 

2.59 

51  .49 

3.94 

9 

20.58 

0.47 

30.68 

1.22 

41.13 

2.13 

51.19 

3.14 

10 

20.51 

0.44 

30.72 

0.80 

41  .23 

1  .86 

50.68 

2.64 

12 

20.48 

0.30 

30.73 

0.65 

40.69 

1.35 

50.97 

1  .90 

'.5 

20.47 

0.22 

30.59 

0.51 

40.88 

0.96 

51.14 

1  .54 

16 

20.52 

0.19 

30.75 

0.46 

40.72 

0.84 

51.38 

1  .  40 

18 

20.49 

0.15 

30.63 

0.38 

40.86 

0.69 

51.00 

1  .04 

20 

20.48 

0.12 

30.66 

0.30 

40.67 

0.61 

51  .02 

1  .00 

a 

20.50 

0.08 

30.70 

0.24 

40.80 

0.44 

51.06 

0.83 

25 

20.49 

0.08 

30.66 

0.21 

40.80 

0.42 

51.14 

0.66 

30 

20.50 

0.06 

30.68 

0.16 

40.78 

0.32 

51.01 

0.48 

32 

20.50 

0.04 

30.64 

0.14 

40.87 

0.24 

51.02 

0.41 

35 

20.49 

0.04 

30.66 

0.13 

40.83 

0.23 

51.04 

0.38 

36 

20.49 

0.04 

30.65 

0.13 

40.83 

0.22 

50.97 

0.34 

«0 

20.49 

0.03 

30.66 

0.09 

40.85 

0.19 

51.01 

0.31 

*5 

20.49 

0.02 

30.67 

0.07 

40.85 

0.1  4 

51  .00 

0.26 

«8 

20.49 

0.02 

30.66 

0.06 

40.85 

0.12 

51  .02 

0.22 

50 

20.49 

0.02 

30.66 

0.06 

40.86 

^.12 

51.06 

0.20 

T-(60, 

12.2) 

T-(80, 

16,2) 

T-O00 

>.20.2) 

T- ( 1 20 

1,24,2) 

Points 

8-61 

.22 

8  -  81. 

61 

8  -  102.00 

R  -  122.39 

i  Traok 

I 

a 

X 

s 

X 

s 

X 

3 

1 

67.16 

22.79 

114.11 

228.2 

139.77 

214.9 

275.42 

825.5 

5 

63.25 

14.24 

90.03 

30.21 

120.45 

53.90 

164.73 

232.6 

6 

61  .76 

8.41 

87.36 

20.08 

112.51 

50.76 

157.78 

71  .48 

8 

61  .85 

7.10 

81.91 

11.84 

104.50 

17.73 

123.10 

24.89 

9 

60.93 

4.61 

83.36 

8.78 

104.73 

18.15 

125.84 

21  .23 

10 

62.08 

4.56 

82.36 

8.29 

102.57 

10.1  1 

125.09 

18.37 

12 

60.72 

2.86 

82.29 

5.85 

103.77 

9.49 

123.08 

12.99 

15 

61  .54 

2.28 

82.00 

4.35 

102.68 

6.02 

121.22 

9.14 

16 

61  .25 

1  .94 

81  .96 

3.82 

102.56 

6.54 

121.17 

8.56 

18 

61 .13 

1  ,70 

81  .55 

2.91 

102.28 

4.95 

123-65 

8.13 

20 

61.41 

1  .53 

81  .59 

2 . 6 1 

102.08 

3.80 

122.39 

5.63 

20 

61.12 

0.98 

81  .95 

2.05 

102.46 

3.14 

122.22 

4.84 

25 

61.11 

1  .08 

81  .64 

1 .89 

102.19 

2.58 

122.62 

4.30 

30 

61.27 

0.80 

81  ,48 

1.43 

102.41 

2.19 

122.25 

3.09 

32 

61  .23 

0.62 

81  .60 

1  .24 

102.29 

1.92 

122.44 

2.91 

35 

61.16 

0.59 

81.82 

1  .01 

102.10 

1  .77 

122.66 

2.46 

36 

61.18 

0.58 

81  .66 

1  .20 

102.11 

1.65 

122.37 

2.35 

*0 

61  .26 

0.45 

81 .67 

0.83 

101.96 

1.18 

122.76 

1  .81 

#5 

61  .23 

0.38 

81  .65 

0.73 

102.26 

1 .21 

122.43 

1.81 

48 

61  .22 

0.34 

81  .62 

0.64 

102.14 

1  .03 

122.33 

1  .57 

50 

61.21 

0.31 

81  .64 

0.56 

102.07 

0.99 

122.38 

1  .34 

61 


table  A 1 0 


Least  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Cell  a  -  0.333  nrad 

Initial  Observer  Veloolty:  <.i6,0,o>  nm/s 
Observer  Acceleration t  <. 006,0, 0>  no/s1 


T-(20. 

*.2) 

T-(30. 

6.2) 

T-(40, 

8,2) 

T  •  (  50  , 

10,2) 

Points 

R_«  20 

.49 

R_-  30 

.66 

R  _•  40 

.84 

R_-  40 

.03 

Track 

X 

s 

X 

s 

X 

3 

X 

3 

4 

20.82 

1  .72 

31  .64 

4.06 

41  .59 

7.26 

52.18 

13.67 

5 

20.  bo 

1  .05 

30.91 

2.77 

42.18 

5.42 

53.22 

9.42 

6 

20.56 

0.85 

30.80 

2.01 

40.84 

4.17 

52.19 

5.94 

8 

20.46 

0.46 

30.65 

1  .29 

40.91 

2.26 

51  .  40 

3.46 

9 

20.57 

0.40 

30.67  • 

’  1  .04 

41  .07 

1  .83 

51.16 

2.73 

10 

20.51 

0.37 

30.70 

0.67 

41.16 

1  .57 

50.69 

2.25 

12 

20.48 

0.25 

30.71 

0.53 

40.71 

1.11 

50.97 

1  .55 

15 

20.47 

0.16 

30.61 

0.39 

40.87 

0.75 

51.11 

1  .20 

16 

20.51 

0.14 

30.73 

0.35 

40.74 

0.64 

51  .28 

1  .07 

18 

20.49 

0.10 

30.64 

0.28 

40.84 

0.50 

51.00 

0.78 

20 

20.49 

0.08 

30.66 

0.21 

40.73 

0.43 

51.01 

0.71 

2t 

20.50 

0.05 

30.68 

0.16 

40.81 

0.29 

51  .05 

0.55 

25 

20.49 

0.05 

30.66 

0.13 

40.81 

0.27 

51.10 

0.43 

30 

20.50 

0.03 

30.67 

0.09 

40.81 

0.19 

51  .01 

0.28 

32 

20.50 

0.02 

30.65 

0.07 

40.85 

0.13 

51.02 

0.24 

35 

20.49 

0.02 

30.66 

0.06 

40.84 

0.12 

51  .04 

0.21 

36 

20.49 

0.02 

30.65 

0.06 

40.84 

o.n 

51  .00 

0.18 

bo 

20.49 

0.01 

30.66 

0.04 

40.84 

0.09 

51  .02 

0.15 

b5 

20.49 

0.01 

30.66 

0.03 

40.84 

0.06 

51.02 

0.12 

b8 

20.49 

0.00 

30.66 

0.02 

40.84 

0.05 

51  .02 

0.09 

50 

20.49 

0.00 

30.66 

0.02 

40.85 

0.05 

51  .04 

0.08 

T-C60, 

12.2) 

T-(80. 

16,2) 

T-(100 

.20.2) 

T • ( 1 20 

.24,2) 

Points 

R_-  61 

.22 

R  -  81  . 

61 

R  -  102.00 

R  -  122.39 

i  Track 

X 

s 

X 

s 

X 

s 

X 

s 

b 

66.41 

20.54 

*  138.64 

419.9 

145.12 

313-4 

317.74 

• 

5 

62.80 

12.78 

88.74 

26.72 

117.43 

44.81 

150.68 

148.8 

6 

61.65 

7.60 

86.63 

17.96 

109.83 

39.41 

150.72 

55.97 

8 

61 .68 

6.08 

81  .63 

10.02 

103.82 

15.25 

122.38 

21  .42 

9 

60.92 

3.96 

83.02 

7.54 

104.04 

15.14 

125.08 

18.06 

10 

61  .92 

3.81 

82.11 

6.88 

102.34 

8.67 

124.37 

15.27 

12 

60.81 

2.37 

82.13 

4.78 

103.28 

7.71 

122.84 

10.66 

15 

61.46 

1 .76 

81  .89 

3.40 

102.49 

4.67 

121.32 

7.15 

16 

61 .24 

1.49 

81  .84 

2.93 

102.39 

4.97 

121.38 

6.58 

18 

61.15 

1  .27 

81  .56 

2.18 

102.17 

3.71 

123-23 

6.00 

20 

61 .36 

1.10 

81  .57 

1  .92 

102.02 

2.73 

122.36 

4.09 

2b 

61.14 

0.66 

81 .82 

1 .38 

102.29 

2.14 

122.23 

3.28 

25 

61.15 

0.72 

81  .61 

1 .27 

102.12 

1  .75 

122.50 

2.82 

30 

61.24 

0.49 

81  .53 

0.88 

102.26 

1  .37 

122.26 

1  .88 

32 

61  .22 

0.37 

81  .60 

0.74 

102.19 

1.17 

122.41 

1  .77 

35 

61.18 

0.32 

81.72 

0.57 

102.04 

1 .00 

122.51 

1.44 

36 

61  .CO 

0.32 

81  .64 

0.67 

102.06 

0.92 

122.36 

1.35 

bo 

61.25 

0.23 

81  .65 

0.45 

101.99 

0.63 

122.59 

0.97 

b5 

61  .23 

0.18 

81  .62 

0.36' 

102.13 

0.58 

122.39 

0.90 

48 

61 .22 

0.15 

81  .62 

0.30 

102.06 

0.49 

122.36 

0.77 

50 

61.22 

0.13 

81  .61 

0.25 

102.02 

0.45 

122.39 

0.65 
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Table  At  1 


Least  Squares  Estlaates  or  Fixed  Target  Position 
100  Saaplea  per  Cell  9  ■  .333  arad 

Initial  Observer  Velooity:  <.l6,0,0> 
Observer  Acceleration:  <0,-.006,0> 


T-(20, 

4,2) 

T-(30. 

6.2) 

T«(40.a 

i.2) 

T-(50, 

10.2) 

Points 

R_-  20 

.49 

R_*  30 

.66 

R_«  40. 

84 

R 40 

.03 

Track 

X 

a 

X 

a 

X 

a 

X 

a 

4 

20.74 

1  .50 

31.38 

3.3“ 

41.37 

5.97 

51  .50 

10.28 

5 

20.40 

0.86 

30.84 

2.21 

41.82 

4.19 

52.36 

7.22 

6 

20.53 

0.67 

30.75 

1  .50 

40.77 

3.09 

51.81 

4.31 

8 

20.46 

0.3“ 

30.65 

0.90 

40.84 

1  .59 

51  .21 

2.39 

9 

20.54 

0.28 

30.65 

0.72 

40.96 

1  .22 

51.15 

1  .85 

1  0 

20.51 

0.24 

30.68 

0.44 

41  .04 

1  .02 

.  50.78 

1  . 46 

12 

20.48 

0.16 

30.70 

0.33 

40.75 

0.68 

51  .00 

0.95 

15 

20.48 

0.10 

30.63 

0.23 

40.84 

0.43 

51  .08 

0.67 

16 

20.50 

0.08 

30.70 

0.20 

40.77 

0.36 

51.15 

0.59 

18 

20.49 

0.06 

30.65 

0.16 

40.85 

0.27 

51  .01 

0.42 

20 

20.49 

0.05 

30.66 

0.11 

40.78 

0.23 

51  .01 

0.36 

24 

20.49 

0.03 

30.67 

0.08 

40.83 

0.14 

51  .04 

0.26 

25 

20.49 

0.03 

30.66 

0.07 

40.83 

0.13 

51.06 

0.20 

30 

20.49 

0.02 

30.67 

0.05 

40.82 

0.09 

51  .02 

0.13 

32 

20.50 

0.01 

30.66 

0.04 

40.85 

0.06 

51  .03 

0.11 

35 

20.49 

0.01 

30.66 

0.03 

40.84 

0.06 

51  .03 

0.09 

36 

20.49 

0.01 

30.66 

0.03 

40.84 

0.05 

51  .01 

0.08 

40 

20.49 

0.01 

30.66 

0.02 

40.84 

0.04 

51  .03 

0.06 

45 

20.49 

0.01 

30.66 

0.02 

40.84 

0.03 

51  .02 

0.05 

48 

20.49 

0.01 

30.66 

0.01 

40.84 

0.02 

51  .03 

0.04 

50 

20.49 

0.00 

30.66 

0.01 

40.84 

0.02 

51  .03 

0.04 

T-(60. 

12,2) 

• 

O 

to 

w 

1 

1- 

16,2) 

H 

• 

O 

O 

,20.2) 

T-( 1 20 

.24.2) 

#  Points  R  -  61.22  R  -  81. 6l  R  -  102.00  R  -  122.39 


on  Track  x  s  xs  xs  xs 

4  64 . 68  15.93  90.41  78.83  130.47  230. 0  214.16  644.3 

5  61.90  9.40  85.28  18.25  111.74  30.19  138.15  119.8 

6  61.30  5.54  84,55  12.28  105.27  22.72  139.72  33-28 


8 

61.39 

4.01 

81  .37 

6.48 

102.79 

10.22 

121 .63 

14.31 

9 

60.98 

2.65 

82.37 

4.92 

102.80 

9.24 

123.70 

11.24 

10 

61.64 

2.35 

81  .80 

4.14 

101 .98 

5.58 

123-12 

9.25 

12 

60.96 

1.42 

81  .86 

2.82 

102.59 

4.50 

122.46 

6.25 

15 

61.34 

0.96 

81  .72 

1  .87 

102.24 

2.53 

121 .69 

3-90 

16 

61  .21 

0.79 

81  .70 

1 .58 

102.16 

2.63 

121.75 

3.50 

18 

61.18 

0.67 

81  .59 

1.11 

102.05 

1.91 

122.73 

3.00 

20 

61.29 

0.55 

81  .57 

0.97 

101 .97 

1.35 

122.34 

1  .99 

24 

61.18 

0.32 

81.70 

0.64 

102.13 

0.98 

122.30 

1  .50 

25 

61.19 

0.34 

81  .59 

0.59 

102.04 

0.80 

122.41 

1  .25 

30 

61  .23 

0.22 

81.58 

0.37 

102.11 

0.60 

122.31 

0.78 

32 

61  .22 

0.16 

81  .61 

0.32 

102.09 

0.50 

122.39 

0.75 

35 

61  .20 

0.14 

81  .65 

0.24 

102.01 

0.41 

122.43 

0.59 

36 

61 .11 

0.14 

81  .62 

0.28 

102.03 

0.38 

122.37 

0.54 

40 

61  .23 

0.10 

81  .63 

0.18 

102.00 

0.25 

122.46 

0.38 

45 

61  .22 

0.08 

81  .61 

0.15 

102.05 

0.23 

122.39 

0.34 

48 

61  .22 

0.07 

81  .61 

0.12 

102.02 

0.19 

122.38 

0.30 

50 

61  .22 

0.06 

81  .61 

0.10 

102.01 

0.17 

122.39 

0 . 24 

63 


Table  A 1 2 


Least  Squares  Eatiaates  of  Fixed  Target  Position 
100  Saaplaa  par  Call  a  -  0.333  arad 

Target  Range  »  to  na  Various  Directions 
Constant  Observer  Velocity*  <.16,0,0>  na/s 


T-(28,3 

.28.3. 

2) 

T-<30. 

6,25. 7, 2> 

« 

V 

i 

6 , 20 . 0 , 2> 

T«<3?.5 

.13.7 

#  Points 

on  Traok 

X 

s 

X 

a 

X 

s 

X 

s 

5 

39  .91 

1  .  t2 

tO. 07 

1 .81 

39.99 

2.26 

tO.  38 

3.10 

8 

to. it 

0.7t 

tO. 08 

0.77 

39.99 

1  .07 

tO. 16 

1  .56 

10 

to. 02 

o.te 

39.99 

0.57 

39.96 

0.76 

to. 06 

1  .07 

12 

39*96 

0.37 

tO.  02 

0 .  t6 

to. 01 

0.52 

39.90 

0.73 

15 

tO. 02 

0.26 

tO. 02 

0.32 

to. 01 

0.39 

39.99 

0.5t 

18 

tO. 02 

0.22 

tO. 05 

0.25 

to. 01 

0.25 

39.99 

0.39 

20 

to. 00 

0.18 

to. 01 

0.17 

tO. 02 

0.21 

39.98 

0.33 

25 

39.99 

0.11 

39.99 

o.it 

tO.  02 

0.17 

to. 00 

0.26 

30 

to. 00 

0.09 

to. 00 

0.10 

39.99 

0.13 

39.98 

0.16 

35 

to. 01 

0.07 

to. 00 

0.07 

39.99 

0.10 

tO. 02 

0.12 

to 

to. 00 

0.05 

to. 00 

0.06 

39.98 

0.08 

39.99 

0.10 

t5 

to. 00 

o.ot 

to. 00 

0.05 

39.99 

0.06 

to. 00 

0.09 

50 

39-99 

o.ot 

to. 00 

o.ot 

'to. 00 

0.05 

to. 00 

0.07 

T-(38.6 

.10.3. 

2) 

T-<39 . 

3.  6.9 ,2> 

T«<39. 

8,  3 . 5 , 2> 

#  Points 

on  Traok 

X 

s 

X 

s 

X 

3 

5 

to. 13 

t.30 

tO. 22 

6 .  t7 

ti  .tt 

12.56 

8 

tO. 05 

1 .98 

t0.5t 

2.97 

tO.  29 

t.to 

10 

tO.  29 

1.28 

39.68 

1  .73 

to. 13 

3.35 

12 

tO. 05 

1  .02 

39.88 

1.29 

tO.  17 

2.53 

15 

39.99 

0.65 

to. 03 

0.91 

tO  .  2t 

1  .8t 

18 

to. 01 

0.58 

to .  1 1 

0.72 

39. 9t 

1  .32 

20 

to. 01 

0.  t2 

39.93 

0.71 

tO. 21 

1  .09 

25 

39.99 

0.30 

39.99 

O.tl 

39.96 

0.77 

30 

39.96 

0.21 

tO. 02 

0.33 

tO. 09 

0.55 

35 

to. 03 

0.18 

39.96 

0.21 

39.93 

0.39 

to 

to. 01 

0.12 

to. 00 

0.21 

to. 01 

0.35 

t5 

39.99 

0.12 

to. 00 

0.17 

39.99 

0.30 

50 

tO.  02 

0.09 

39.98 

0.12 

tO. 02 

0.23 

Tabla  A1 J 

Laaat  Squares  Katiaatas  of  Flxad  Target  Position 
100  Saaplaa  par  Call  «  -  0.333  arad 

Target  Range  -  AO  na  Various  Directions 
Initial  Observer  Velocity!  <.16,0,0>  na/s 
Observer  Acceleration t  <0,-.006,0>  na/a* 


T-(28. 

3.28.3. 

2) 

T-<30. 

6 , 25 • 7  1 2> 

T-<3A.6 

■  20 . 0  i  2> 

T-<37 .5 

.13.7 

#  Points 

on  Track 

X 

a 

X 

a 

X 

a 

X 

a 

5 

AO  28 

1.31 

39.98 

1  .  A 1 

AO.  A5 

1  .56 

AO. It 

2.28 

8 

AO. 18 

0.53 

AO.OA 

0.69 

39.93 

0.7A 

39.92 

0.96 

10 

39. 9A 

0.  AO 

39.89 

0.  A? 

39.97 

0.51 

39.96 

0.75 

12 

AO. 05 

0.28 

AO. 00 

0.29 

AO. 07 

0.36 

AO. 00 

0.  A8 

15 

AO. 05 

0.21 

39.99 

0.21 

AO. 00 

0.23 

AO. 00 

0.31 

18 

39.98 

0.1  A 

AO. 00 

0.16 

39.99 

0.19 

39.99 

0.19 

20 

AO. 01 

0.12 

39.97 

0.13 

39.98 

0.1  A 

39.98 

0.17 

25 

AO. 01 

0.09 

39.99 

0.09 

39.99 

0.10 

AO. 02 

0.11 

30 

AO. 01 

0.05 

39.99 

0.06 

AO. 00 

0.06 

AO. 01 

0.06 

35 

AO. 00 

0.05 

39.99 

O.OA 

AO. 01 

0.05 

AO. 00 

0.05 

AO 

AO.  00 

0.03 

AO. 00 

O.OA 

AO. 00 

O.OA 

AO. 00 

0.03 

A5 

AO. 00 

0.03 

39.99 

0.03 

AO. 00 

0.03 

AO. 00 

0.03 

50 

AO. 00 

0.02 

AO. 00 

0.02 

AO. 00 

0.02 

AO. 00 

0.02 

T-(38. 

6.10.3. 

2) 

T-<39. 

3.  6 .9  »2> 

T-<39 .8 

.  3 • 5 , 2> 

#  Points 

on  Track 

X 

a 

X 

a 

X 

a 

5 

AO .  1  A 

3.10 

A1  .OA 

A. 88 

AO.  10 

6  .OA 

8 

39.79 

1.3A 

AO. 12 

1 .61 

AO. 59 

2.37 

10 

39.93 

0.80 

39.91 

0.98 

AO. 10 

1 ,2A 

12 

39.95 

0.52 

39.99 

0.73 

39.92 

0.87 

15 

39.96 

0.37 

AO. 00 

0.A3 

AO. 03 

0.51 

18 

39.98 

0.25 

39.99 

0.26 

AO. 02 

0.36 

20 

AO. 01 

0.19 

AO. 00 

0.41 

39.99 

0 . 2  A 

25 

AO. 01 

0.12 

39.99 

0 . 1  A 

AO. 02  ~ 

0.13 

30 

AO.  00 

0.08 

AO. 01 

0.07 

AO. 00 

0.09 

35 

39.99 

0.05 

39.99 

0.06 

AO. 00 

0.06 

AO 

AO. 00 

0.05 

AO. 00 

0.03 

AO. 00 

O.OA 

A5 

AO. 00 

0.03 

39.99 

0.03 

AO. 00 

0.03 

50 

AO.  00 

0.02 

AO. 00 

0.02 

AO. 00 

0.03 
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Table  A14 


Least  Squares  Estimates  of  Fixed  Target  Position 
100  Samples  per  Cell  o  -  0.333  mrad 

far  get  Location:  (  50,40,2)  Range  -  64.06  nm 

'Initial  Observer  Velocity:  <.16,0,0>  nm/s 


Poi nts 

a«< .006 

,0,0> 

?"<0 

,0,0> 

a»<0 , 

.006 ,0> 

a- <0 ,  - 

.006 ,0> 

l  Track 

X 

s 

X 

s 

X 

3 

X 

3 

4 

64.32 

5.58 

64.74 

6.31 

63.12 

6.87 

64.18 

5  .58 

5 

64.82 

4.24 

64.11 

4.18 

64.48 

4.71 

63.74 

4 .03 

6 

63.97 

3.26 

63.95 

2.83 

64,63 

3.83 

65.28 

2.93 

8 

64.05 

1  .87 

64.09 

2.28 

63.86 

2.43 

63.81 

1  .80 

9 

64.22 

1.48 

63.91 

1  .65 

64.33 

2.01 

64.12 

1  .39 

10 

64.29 

1  .29 

64.30 

1  .50 

64.15 

1  .87 

64.10 

1  .20 

12 

63.95 

0.93 

63.86 

1  .01 

64.15 

1  .51 

64.04 

0 , 86 

15 

64.06 

0.64 

64.16 

0.79 

64.12 

1  .21 

63.95 

0. 60 

16 

63.96 

0.54 

64.06 

0.66 

64.13 

1.12 

63.95 

0.54 

18 

64.08 

0,44 

64.02 

0.59 

64.01 

0.88 

64.11 

0 . 49 

20 

63.95 

0.38 

64.12 

0.54 

64.05 

0.85 

64.05 

0.33 

24 

64.04 

0.26 

64.03 

0.35 

64.19 

0.78 

64.04 

0 . 27 

25 

64.04 

0.24 

64.02 

0.39 

64.09 

0.75 

64.07 

0.23 

30 

64.03 

0.17 

64.08 

0.29 

63.98 

0.72 

64.05 

0.15 

32 

64.07 

0.13 

64.07 

0.22 

64.05 

0.69 

64.06 

0.15 

35 

64.06 

0.11 

64.04 

0.22 

64.18 

0.66 

64.07 

0.12 

36 

64.06 

0.11 

64.05 

0.21 

64.08 

0.69 

64.06 

0.11 

40 

64.07 

0.09 

64.08 

0.16 

63.95 

0.60 

64.08 

0 .08 

45 

64.06 

0.06 

64.06 

0.14 

64.10 

0.52 

64.06 

0.08 

48 

64.07 

0.05 

64.06 

0.13 

64.01 

0.40 

64.06 

0.07 

50 

64.07 

0.05 

64.06 

0.11 

64.08 

0.30 

64.06 

0.06 
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Table  A15 


Least  Squares  Estimates  of  Fixed  Target  Position 


100  Samples 

per  Cell 

o  ■ 

0.333  mrad 

Target  Location:  (50 

,10,2) 

Range 

-  51.03  nm 

Initial 

Observer 

Velocity: 

< .  1 6 , 0 

,0>  nm/s 

a«<.006,0,0> 

a«<0 

,0,0> 

a-<0 , 

.006 ,0> 

a-<0,- 

.006,0> 

Points 

i  Track 

X 

s 

X 

s 

X 

s 

X 

s 

4 

56.52 

23.67 

55.57 

16.03 

56.92 

35.39 

51 .27 

10.87 

5 

51.20 

8.09 

52.53 

8.86 

57.59 

18.41 

52.46 

6.11 

6 

51.77 

6.48 

51 .90 

6.99 

54.29 

12.98 

51.30 

4.89 

a 

50.82 

3.42 

51 .11 

4.27 

53.92 

12.09 

51.10 

2.15 

9 

51 .68 

2.93 

51 .14 

3.69 

50.85 

8.21 

51 .13 

2.15 

10 

51.30 

2.56 

51 .24 

2.40 

52.64 

10.01 

50.99 

1.37 

12 

50.95 

1 .82 

51.30 

1 .93 

51 .19 

5.91 

51.15 

1  .08 

15 

50.88 

1.26 

50.87 

1.52 

50.86 

3.21 

51  .08 

0.61 

16 

51.21 

1 .09 

51 .35 

1 .39 

51.03 

2.38 

51 .06 

0.62 

18 

51 .04 

0.82 

50.96 

1.16 

51 .21 

1  .59 

51.04 

0.46 

20 

50.96 

0.66 

51.05 

0.95 

51.16 

1  .06 

51 .02 

0.32 

24 

51.03 

0.43 

51.17 

0.72 

51 .01 

0.61 

51 .06 

0.23 

25 

51 .00 

0.45 

51 .05 

0.65 

50.99 

0.46 

51.04 

0.19 

30 

51.04 

0.32 

51 .08 

0.52 

51 .05 

0.24 

51 .05 

0.14 

32 

51.06 

0.23 

50.98 

0.44 

51 .04 

0.20 

51 .05 

0.12 

35 

51.03 

0.22 

51 .04 

0.41 

51 .02 

0.16 

f  .03 

0.09 

36 

51.02 

0.19 

51.00 

0.42 

51 .05 

0.14 

51 .03 

0.09 

40 

51.02 

0.15 

51 .04 

0.29 

51.03 

0.10 

51.03 

0.06 

45 

51.02 

0.10 

51.07 

0.25 

51.04 

0.08 

51.04 

0.05 

48 

51.0? 

0.10 

51 .04 

0.20 

51.03 

0.06 

51.03 

0.04 

50 

51.02 

0.07 

51 .02 

0.22 

51 .02 

0.05 

51.03 

0.04 
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TABLE  A.  1  6 


Azimuth  Only  Least  Squares  Estimates  of  Fixed  Target  Position 


100 

Samples 

per  Cell 

a  * 

0.333  mr ad 

Constant 

Observer 

Velocity: 

:  <  .  1 6 , 0 , 0>  nm/s 

T-(20, 

16,2) 

T-(30, 

24,2) 

T-(40.32.2) 

T-(50,40,2) 

Points 

R_-  25 

.69 

R_«  38 

.47 

R_-  51 

.26 

R_-  64 

.06 

Track 

X 

s 

X 

s 

X 

s 

X 

3 

4 

25.70 

1  .00 

38.55 

2.15 

50.69 

3.58 

62.59 

5.47 

5 

25.59 

0.65 

38.54 

1.54 

51.35 

2.92 

64.06 

4.69 

6 

25.65 

0.53 

38.39 

1  .09 

50.89 

2.24 

63.85 

3.02 

8 

25.63 

0.31 

38.36 

0.76 

51.06 

1  .38 

63.73 

2.02 

9 

25.72 

0.26 

38.38 

0.66 

51.21 

1  .08 

63.85 

1  .63 

1  0 

25.70 

0.23 

38.44 

0 .43 

51.30 

0.96 

63.58 

1  .28 

1  2 

25.65 

0.17 

38.48 

0.35 

51  .09 

0.71 

63.87 

0.99 

15 

25.67 

0.12 

38.42 

0.27 

51  .20 

0.52 

63-99 

0 . 80 

16 

25.70 

0.11 

38.51 

0.25 

51.13 

0 . 45 

64.12 

0.73 

18 

25.68 

0.08 

38.44 

0.22 

51.24 

0 . 38 

63.95 

0.54 

20 

25.68 

0.07 

38.46 

0.18 

51.13 

0.33 

63.96 

0.52 

24 

25.69 

0.05 

38.49 

0.13 

51  .22 

0.24 

64.02 

0.43 

25 

25.68 

0.05 

38.46 

0.12 

51.22 

0.22 

64.07 

0.35 

30 

25.69 

0.04 

38.47 

0.10 

51.21 

0.18 

64.03 

0.26 

32 

25.69 

0.03 

38.46 

0.08 

51.26 

0.13 

64.03 

0.21 

35 

25.69 

0.03 

38.47 

0.07 

51.25 

0.12 

64.05 

0.21 

36 

25.69 

0.03 

38.46 

0.08 

51.26 

0.12 

64.01 

0.19 

40 

25.69 

0.02 

38.47 

0.05 

51.26 

0.11 

64.04 

0.17 

45 

25.69 

0.02 

38.48 

0.05 

51.26 

0.08 

64.03 

0.15 

48 

25.69 

0.02 

38.47 

0.04 

51  .27 

0.07 

64.05 

0.12 

50 

25.69 

0.01 

38.47 

0.04 

51.27 

0.07 

64.07 

0.11 

T- ( 60 , 

48,2) 

T- ( 80 , 

64,2) 

T-M00 

.80.2) 

T- ( 1 20 

.96.2) 

Points 

R_-  76 

.86 

R  -  102 

.47 

R  -  128.08 

R  -  153.69 

i  Track 

X 

s 

X 

s 

X 

3 

X 

3 

4 

76.51 

8.70 

96.90 

15.80 

120.07 

21  .27 

145.50 

37.44 

5 

75.49 

5.70  • 

100.89 

10.89 

125.27 

16.16 

143.75 

23.39 

6 

75.76 

4.15 

101.71 

8.42 

124.04 

13-60 

153.66 

19.79 

8 

76.24 

3.21 

100.42 

5.10 

125.58 

8.04 

147.52 

11.36 

9 

76.11 

2.40 

101.98 

4.33 

126.07 

7.98 

150.42 

9.66 

10 

76.83 

2.16 

101 .79 

3.73 

126.05 

5.06 

150.46 

8.79 

12 

76.27 

1.43 

102.00 

2.90 

127.30 

4.50 

151 .05 

6.18 

1  5 

76.80 

1.16 

102.10 

2.18 

127.34 

3-03 

1 51 .22 

4.47 

16 

76.66 

0.95 

102,11 

1.95 

127.39 

3.21 

151 .29 

4.20 

18  ' 

76.66 

0.85 

102.08 

1  .41 

127.49 

2.41 

152.98 

4.04 

20 

76.84 

0.78 

102.1 3 

1.33 

127.45 

1  .95 

152.62 

2.80 

24 

76.73 

0.51 

102.42 

1  .03 

127.93 

1  .58 

152.88 

2.44 

25 

76.73 

0.57 

102.26 

0.98 

127.78 

1  .28 

153.07 

2.24 

30 

76.84 

0.42 

102.29 

0.72 

128.05 

1.12 

153.15 

1  .60 

32 

76.83 

0.32 

102.35 

0.64 

128.01 

0.98 

153.27 

1.51 

35 

76.79 

0.32 

102.48 

0.53 

127.95 

0.90 

153.48 

1  .26 

36 

76.80 

0.31 

102.39 

0.63 

127.95 

0.86 

153.35 

1  .  20 

40 

76.85 

0.24 

102.42 

0.42 

127.90 

0.63 

153.59 

0.92 

**5 

76.84 

0.21 

102.43 

0.39 

1 28 . 1  0 

0.63 

153.51 

0.92 

48 

76.84 

0.19 

102.43 

0.33 

128.05 

0.53 

153.49 

0.81 

50 

76.84 

0.17 

102,44 

0.29 

128.03 

0.52 

153-52 

0.64 

TABLE  A  1 7 

Median  Based  Estimates  of  Fixed  Target  Position 
Medians  of  all  Pairs  of  Two-point  Estimates 
100  Samples  per  Cell  o  -  0.333  mrad 


Constant 

Observer 

Velocity : 

<  .  1  6 , 0 

,0> 

T  -  (50,10,2) 

T  -  (50,40, 

#  Points 

R  * 

51.03 

R  ■ 

64.06 

on  Track 

X 

s 

X 

s 

4 

51  .53 

1  5.55 

65.09 

6.30 

5 

49.69 

10.53 

64.23 

4.46 

6 

50.03 

7.99 

64.20 

3.59 

8 

50.38 

4.32 

63.99 

2.29 

9 

50.11 

3.82 

64.02 

1 .92 

10 

50.93 

3.10 

64.15 

1  .38 

1  2 

50.49 

2.31 

64.24 

0.95 

1  5 

50.52 

1  .56 

64.01 

0.85 

1  6 

51.19 

1  .43 

64.20 

0.75 

18 

50.84 

1.12 

64.02 

0.63 

20 

50.82 

0.96 

64.05 

0.53 

24 

50.93 

0.69 

64.13 

0.38 

25 

50.91 

0.74 

64.06 

0.36 

30 

50.97 

0.54 

64.09 

0.28 

32 

51  .07 

0.33 

64.04 

0.25 

35 

51  .01 

0.42 

64.07 

0.23 

36 

50.97 

0.36 

64.05 

0.24 

40 

50.99 

0.32 

64.06 

0.15 

45 

50.98 

0.25 

64.09 

0.14 

48 

51  .01 

0.26 

6  4.07 

0.11 

50 

51  .02 

0.20 

64.06 

0.12 

a 
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TABLE  A  1 8 

Least  Squares  Estimates  of  Constant  Velocity  Target  Parameters 
100  Samples  per  Cell  a  -  0.333  mrad 

Target:  Initial  Position  -  (50,10,1)  Range  -  51.00 

Various  Velocities 

Initial  Observer  Velocity:  <.16,0,0> 

Various  Observer  Accelerations 


a  -  <0 , 

-.002, 

.006> 

a  *  <  0 

, - . 006 , . 

002> 

V  -  <-. 

1  .-.1  , 

0> 

V  -  <- 

. 1 ,-. 1 ,0> 

Range 

Vx 

Range 

Vx 

#  Points 

on  Track 

X 

s 

X 

s 

X 

s 

X 

3 

15 

50.77 

5.37 

-0.05 

0.26 

50.45 

5.05 

0.12 

0.55 

18 

50.25 

4.39 

-0.07 

0.18 

49.03 

4.68 

0.10 

0.35 

20 

50.66 

3.10 

-0.08 

0.12 

49.02 

5.31 

0.06 

0.36 

24 

50.94 

1  .61 

-0.09 

0.08 

50.35 

3.28 

-0.06 

0.16 

30 

50.40 

2.59 

-0.08 

0.07 

50.36 

2.83 

-0.08 

0.09 

36 

50.68 

1  .55 

-0.09 

0.03 

50.77 

1  .93 

-0.10 

0.05 

40 

50.97 

0.95 

-0.10 

0.02 

50.80 

1  .64 

-0.10 

0.04 

48 

50.98 

0.64 

-0.10 

0.01 

51  .03 

0.72 

-0.10 

0.01 

54 

50.97 

0.53 

-0.10 

0.01 

50.81 

1  .28 

-0.10 

0.02 

60 

50.96 

0.34 

-0.10 

0.00 

51  .08 

0.45 

-0.10 

0.01 

75 

50.91 

0.17 

-0.10 

0.00 

51  .01 

0.26 

-0.10 

0.00 

a  •  <0 

,-.006 

, .  002  > 

a  - 

<0 ,-.002 

, . 006> 

V  -  <0 

,  — .  1  6  , 1 

0> 

V  - 

<0 ,  1 6  , 

0> 

Range 

vy 

Range 

Vy 

#  Points 

on  Track 

X 

s 

X 

s 

X 

s 

X 

3 

15 

53.92 

9.81 

-0.14 

0.05 

51  .38 

4.31 

-0.17 

0.02 

18 

50.71 

3.80 

-0.16 

0.02 

51  .06 

2.66 

-0.16 

0.01 

20 

51  .08 

2.88 

-0.16 

0.02 

51  .22 

2.22 

-0.16 

0.01 

24 

51.16 

1  .27 

-0.16 

0.01 

50.74 

1  .40 

-0.16 

0.00 

30 

50.96 

0.82 

-0.16 

0.01 

50.94 

0.87 

-0.16 

0.00 

36 

51.16 

0.58 

-0.16 

0.00 

50.84 

0.64 

-0.16 

0.00 

40 

50.96 

0.49 

-0.16 

0.00 

50.93 

0.45 

-0.16 

0.00 

48 

50.99 

0.45 

-0.16 

0.00 

51  .05 

0.30 

-0.16 

0.00 

54 

50.87 

0.35 

-0.16 

0.00 

50.98 

0.28 

-0.16 

0.00 

60 

50.88 

0.30 

-0.16 

0.00 

50.96 

0.22 

-0.16 

0.00 

75 

50.99 

0.19 

-0.16 

0.00 

50.95 

0.14 

-0.16 

0.00 
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T  A  E  L  E  A  1 '? 

Least  Sauares  Estimates  cf  Constant  Velocity  Target  Parameters 
IOC  Samples  per  Cell  *  0.333  mrac 

Target:  Initial  Position  =  (50,40,1)  Aanqe  =  <*  4 .  C  4 

Velocity :<-.1,-.  1  ,C> 

Initial  Observer  Velocity:  <.16,0,0> 

Various  Ocserver  Accelerations 


a  =  <  G* 

•  006,  . 

OOO 

a  =  <C 

,— .C06,. 

C  Ca> 

Par  ge 

Vx 

flange 

Vx 

Points 

i  Track 

7 

s 

X 

s 

X 

s 

7 

s 

15 

65.33 

5.61 

-C.16 

0.17 

66.23 

8.80 

-C.16 

C.1  9 

18 

65.08 

3.72 

-0.1  5 

C.10 

66.97 

7.33 

-C.15 

0.14 

20 

64.80 

3.C5 

-0.14 

0.08 

64.91 

5.18 

-C.12 

C.C9 

24 

64.14 

2.24 

-0.1C 

C  .  Q  6 

64. 7C 

3.79 

-C.11 

C.C6 

30 

64.47 

1.55 

-0.11 

0.04 

64.3  6 

2.51 

-C.11 

0.C4 

36 

63.  $9 

1.15 

-0.10 

C  .02 

63.86 

1.75 

-C.1C 

G.C2 

40 

63.91 

1.11 

-0.1C 

C  .02 

6C.25 

4.38 

-C.C5 

0.C5 

48 

63.99 

0.7C 

-0.10 

G.01 

63.35 

1.01 

-C.10 

O.C1 

54 

64.09 

C.  64 

-0.1C 

C.00 

64. 1C 

o.ao 

-C.10 

0.C1 

6C 

63.93 

C.44 

-0.10 

c.oo 

64.07 

C.',7 

-C.10 

0.C1 

75 

63.96 

0.20 

-0.1C 

C.00 

63.88 

0.37 

-C.10 

o.:c 

a  s  <0,-002, 

•  0G6> 

a  = 

<0,.00ft. 

.  C02> 

# 

.large 

Vx 

flange 

Vx 

Points 

Track 

K 

s 

7 

s 

X 

s 

7 

s 

15 

64.09 

7.C6 

o.oc 

C.09 

65.80 

8.67 

-C.02 

O.C5 

18 

63.11 

5. 1C 

-0.03 

C.09 

65.68 

6.03 

-C.02 

C.C6 

20 

62.74 

4.24 

-0.C5 

C.08 

65.09 

4.47 

-C.01 

C.C3 

24 

64.39 

2.98 

-0.11 

C.  06 

63. 7C 

2.60 

-C.02 

C.C6 

30 

64.02 

1.92 

-0.10 

C.03 

62.32 

1.72 

-C.01 

O.C5 

36 

64.19 

1.39 

-0.1C 

0.02 

60.74 

0.66 

C.02 

0.C1 

4Q 

64.18 

1.01 

-0.1C 

G.01 

61.42 

2.06 

-C.03 

0.C6 

48 

64.04 

0.79 

-0.10 

C.01 

63.76 

1.00 

-C.10 

0.C2 

54 

63.81 

1.34 

-0.1C 

C.C2 

61.72 

3.34 

-C.06 

G.C5 

60 

64.01 

0.44 

-0.1C 

G.OO 

64.08 

0.52 

-C.10 

C.C1 

75 

63.85 

1.25 

-0.1C 

C.01 

64.02 

0.29 

-C.10 

O.CC 
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TAELe  A1«  (Continued) 


Least  Squares  Estimates  cf  Constant  Velocity  Target  Parameters 
IOC  Saircles  oer  Cell  =  0.333  (r.rac 

Target:  Initial  Position  =  (50/40*1)  Range  =  e  4 . 04 

Velocity:<-.1/-.1  *C> 

Initial  Ctserver  Velocity:  <.16/0/0 
Various  Observer  Accelerations 


a  —  <  0  / 

-.002/  . 

QC6> 

a  =  <3 

/- .  GC6  / 

.  00  2  > 

Range 

Vx 

Range 

Vx 

Points 

i  Track 

7 

s 

7 

s 

7 

s 

X 

s 

15 

67. C7 

9.70 

-C.18 

C.22 

70.  S3 

28.36 

-C.  ’C 

c.45 

18 

64.75 

5.31 

-0.11 

C.11 

64.74 

10. *4 

-C.11 

0.16 

20 

65.32 

5.49 

-0.1? 

0.10 

63.04 

1  Q.  2  A 

-C.08 

Q.16 

24 

64.31 

3.C7 

-0.11 

C.05 

58.43 

6.47 

-C  .02 

O.CR 

30 

63.98 

2.13 

-0.10 

C.  03 

56.56 

1.83 

-C.01 

0.C2 

36 

64.17 

1  .52 

-0.10 

C  .02 

55.47 

1.59 

C.OO 

0.  C  2 

40 

63.98 

1.35 

-0.1C 

C.C2 

55.21 

3.25 

C.CO 

3.C3 

48 

64.  C6 

0.S2 

-0.1C 

C.01 

60. 5C 

5.69 

-C.07 

u.C6 

54 

64. C4 

0.62 

-0.1C 

C.01 

61.86 

4.45 

-C  .08 

0.0  4 

60 

63.93 

0.65 

-0.1C 

C.01 

64.02 

1.04 

-C.10 

0.01 

75 

63.97 

0.21 

-0.1C 

c.oo 

64.  OC 

0.62 

-C.10 

c.cc 

a  *  <Q 

/  a  004  /  • 

QC4> 

Range 

Vx 

Points 

Track 

7 

s 

X 

s 

15 

65.16 

3.7C 

-0.16 

C.  1  5 

18 

64.  32 

5.54 

-0 .16 

C.11 

20 

64.91 

4  •  C  4 

-0.12 

C.O? 

24 

64. C5 

2.54 

-0.11 

C  .06 

30 

63. 85 

1.76 

-0.1C 

0.04 

36 

64.  CO 

1.25 

-0.1C 

C.C2 

40 

64. 04 

0.98 

-0.1C 

0.02 

43 

62.14 

3.23 

-0.07 

C.05 

54 

64.  C5 

0.64 

-0.1C 

C.01 

60 

64.08 

0.52 

-0.1C 

C.01 

75 

63.55 

0.21 

-0.1C 

C.OO 

i 
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TA6LE  A19  (Continued) 


Least  Squares  Estimates  cf  Constant  Velocity  Target  Parameters 
IOC  Samples  cer  Cell  *  0.233  mrac 

Target:  Initial  Position  =  (5C*4Q#1)  Range  =  64.04 

Velocity: <C*-. 16  *Q> 

Initial  Observer  Velocity:  <.16/Q#G> 

Various  Observer  Accelerations 


a  =  <0/ 

* . 006  0 . 

0C2> 

3  =  <C, 

-.CC2,. 

CC6> 

Rarge 

v  y 

Range 

Vy 

Points 

Track 

“x 

s 

7 

s 

7 

s 

X 

s 

15 

70.62 

25.56 

-0.41 

1.33 

73.84 

15.63 

-C.41 

0.4  2 

13 

6  6.85 

6.95 

-0.26 

0.1  R 

67.57 

7.98 

-C.26 

C.2.2 

20 

66.36 

6  .  C  0 

-0.23 

C.15 

63.16 

8.98 

-C.25 

C.2C 

24 

64.72 

2.47 

-0.16 

C.Q4 

64.55 

5. S3 

-G.17 

0.12 

3G 

62.60 

2.24 

-0.14 

C  •  04 

63.95 

3.41 

-C.16 

C.C5 

30 

61.26 

2.5C 

-C.12 

0.04 

63.76 

1.36 

-C.16 

Q.C2 

40 

61.68 

2.65 

-0.13 

0.03 

63.02 

1.35 

-C.15 

Q .  C  2 

43 

61.83 

3.77 

-0.14 

C.04 

63.18 

1.35 

-C.15 

0.C1 

54 

64.90 

0.56 

-0.17 

0.01 

64.94 

0.94 

-C.16 

0.C1 

60 

63.16 

0.82 

-0.15 

0.01 

63.38 

0.51 

-C.16 

o.cc 

75 

63.29 

2.61 

-0.16 

C.02 

63.73 

0.38 

-C.16 

o.cc 
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TABLE  A19  (Continued) 

Least  Squares  Estimates  of  Constant  Velocity  Target  Parameters 
100  Samples  per  Cell  a  -  0.333  mrad 

Target:  Initial  Position  -  (50,40,1)  Range  -  64.04 

Velocity 2 2  , 0> 

Initial  Observer  Velocity:  <.16,Q,0> 

Various  Observer  Accelerations 


a  «  <0 , 

.006 , .006> 

a  -  <0 

, - . 006 ,  . 

006  > 

Range 

Vx 

Range 

VX 

#  Points 

on  Track 

X 

s 

X 

s 

X 

s 

X 

s 

15 

65.42 

6.09 

-0.26 

0.23 

65.57 

9.88 

-0.25 

0.25 

18 

64.44 

3.98 

-0.20 

0.1  6 

65.81 

7.14 

-0.24 

0.15 

20 

63.95 

4.51 

-0.22 

0.14 

63.80 

5.61 

-0.20 

0.12 

24 

64.45 

2.90 

-0.21 

0.08 

64.56 

4.54 

-0.21 

0.08 

30 

64.34 

2.03 

-0.21 

0.04 

64.07 

2.82 

-0.20 

0.04 

36 

64.20 

nJ  .46 

-0.20 

0.02 

63.94 

1  .76 

-0.20 

0.02 

40 

63.97 

1.10 

-0.20 

0.02 

61  .24 

6.01 

-0.16 

0.08 

48 

64.03 

0.70 

-0.20 

0.01 

62.91 

5.55 

-0.19 

0.07 

54 

64.18 

0.59 

-0.20 

0.01 

63.92 

0.73 

-0.20 

0.01 

60 

64.03 

0.39 

-0.20 

0.00 

64.14 

0.52 

-0.20 

0.01 

75 

63.94 

0.17 

-0.20 

0.00 

64.02 

0 . 32 

-0.20 

0.00 

a  •  <0 

,  .002 ,  . 

006> 

a  ■ 

<0 , .006  , 

.  002> 

Range 

Vx 

Range 

Vx 

#  Points 

on  Track 

X 

s 

X 

s 

X 

3 

X 

3 

15 

59.88 

7.65 

0.03 

0.03 

64.38 

8.25 

-0.03 

0.06 

18 

59.86 

4.20 

-0.05 

0.13 

65.52 

6.47 

0.00 

0.02 

20 

60.61 

4.69 

-0.07 

0.13 

63.68 

3.89 

-0.03 

0.07 

24 

63.12 

3.56 

-0.17 

0.08 

61  .91 

2.55 

-0.04 

0.10 

30 

61.37 

4.73 

-0.14 

0.10 

59.81 

2.39 

-0.02 

0.08 

36 

64.02 

1  .43 

-0.20 

0.02 

56.13 

2.20 

0.02 

0.06 

40 

63.87 

1.18 

-0.20 

0.02 

61  .07 

4.45 

-0.13 

0.10 

48 

64.12 

0.78 

-0.20 

0.01 

62.28 

4.42 

-0.17 

0.07 

54 

64.08 

0.51 

-0.20 

0.01 

64.04 

0.76 

-0.20 

0.01 

60 

64.00 

0.48 

-0.20 

0.00 

64.13 

0.58 

-0.20 

0.01 

75 

63.95 

0.28 

-0.20 

0.00 

64.07 

0.32 

-0.20 

0.00 
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TABLE  A20 


Least  Squares  Estimates  of  Constant  Velocity  Target  Position 
tOOO  Samples  per  Cell  o  -  0.333  mr  40  Points/Track  10  Replicates 
Target  —  Location  :  (50 , 40  ,'J)  Range:  64.04  Velocity:  <- .  20  ,  - .  20 , 0> 
Observer  —  Initial  Veloel ty : < . 1 6 , 0 , 0>  Accelerat Ion : <0 , - . 006 , . 006 > 


R 

61  .76 

62.00 

61.88 

61  .87 

61  .  73 

61.50 

61.87 

61 .73 

61  .69 

6  1.9 

s 

5.69 

5.56 

5.64 

5.76 

5.82 

5.96 

5.71 

5.77 

5.33 

5 . 7, 

mln 

48.43 

48.46 

48.54 

48.50 

48.39 

48.51 

48.42 

48.56 

48.47 

4  H  i: 

15 

48.66 

48 ;  67 

48.69 

48.64 

*48.70 

48.68 

48.73 

48.70 

48.67 

“5.7 

55 

48.92 

48.96 

48.94 

48.90 

48.92 

48.92 

48.95 

48.88 

48.92 

-5.9 

105 

49.12 

49.16 

49.08 

49 .09 

49.09 

49.08 

49.12 

49.08 

49.08 

-  9  . ' 

155 

49.41 

59.60 

49.49 

49.41 

49.34 

49.26 

49.42 

49.34 

49.35 

205 

62.03 

62.25 

62.16 

62.17 

62.08 

61  .53 

62.13 

61.80 

61.78 

62.0 

255 

62.61 

62.70 

62.63 

62.64 

62.69 

62.48 

62.64 

62.53 

62.45 

62.6 

305 

62.94 

63.04 

62.95 

62.93 

63.00 

62.78 

62.99 

63.02 

62.77 

63.0 

355 

63.20 

63.26 

63.26 

63.25 

63.30 

*63-09 

63.28 

63.28 

63-17 

63.3 

405 

63.41 

63.48 

63.53 

63.51 

63.45 

63.35 

63.48 

63.45 

63-«9 

63.5 

455 

63.65 

63.73 

63.74 

63-76 

63.69 

63-59 

63.68 

63.68 

63.70 

63.7 

505 

63.85 

63.94 

63.92 

63-93 

63.88 

63.79 

63.92 

63.85 

63.88 

63.9 

555 

64.03 

64.11 

64.10 

64.07 

64.05 

63.95 

64.09 

63.99 

64.04 

64.0 

605 

-  64.18 

64.27 

64.28 

64.26 

64..  29 

64.15- 

64.27 

64.15 

64.24 

6  4.2 

655 

64.37 

64.47 

64.44 

64.45 

64.47 

64.36 

64. 52 

64.42 

64.44 

64.9 

70S 

64.54 

64.71 

64.64 

64.79 

64.66 

64.60 

64.77 

64.59 

64.71 

6  -  .  1 

755 

64.89 

64.95 

64.92 

65.05 

64.89 

64.87 

65.02 

64.88 

64.93 

6  4  .  t 

805 

65.12 

65.20 

65.17- 

65.34 

65.19 

65.17 

65.24 

65.24 

65.21 

65.2 

855 

65.46 

65.48 

65.50 

65.66 

65.49 

65.45 

65.57 

65.57 

65.52 

65.4 

905 

65.74 

65.90 

65.86 

66.04 

65.87 

65.87 

65-95 

65.92 

65.93 

65.9 

955 

66.25 

66.44 

66.33 

66.59 

66.39 

66.48 

66.49 

66.36 

66.1.3 

66.5 

995 

67.39 

67.52 

67.26 

67.57 

67.45 

67.59 

67.53 

67.49 

67.4fc 

6-. 4 

max 

68.78 

68.56 

68.73 

69.22 

69.75 

68.50 

69.26 

68.35 

68.62 

*  G  1 

v7  »  ■ 

* 


75 


0>  Oi  <*> 


